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Abstract 

The effort to resolve data quality issues and improve on the initial data evaluation methodologies of the Sea- 
viewing Wide Field-of-view Sensor (SeaWiFS) Project was an extensive one. These evaluations have resulted, 
to date, in three major reprocessings of the entire data set where each reprocessing addressed the data quality 
issues that could be identified up to the time of each reprocessing. The number of chapters (21) needed to 
document this extensive work in the SeaWiFS Postlaunch Technical Report Series requires three volumes. The 
chapters in Volumes 9, 10, and 11 are in a logical order sequencing through sensor calibration, atmospheric 
correction, masks and flags, product evaluations, and bio-optical algorithms. The first chapter of Volume 9 is 
an overview of the calibration and validation program, including a table of activities from the inception of the 
SeaWiFS Project. Chapter 2 describes the fine adjustments of sensor detector knee radiances, i.e., radiance 
levels where three of the four detectors in each SeaWiFS band saturate. Chapters 3 and 4 describe the analyses 
of the lunar and solar calibration time series, respectively, which are used to track the temporal changes in 
radiometric sensitivity in each band. Chapter 5 outlines the procedure used to adjust band 7 relative to band 8 
to derive reasonable aerosol radiances in band 7 as compared to those in band 8 in the vicinity of Lanai, Hawaii, 
the vicarious calibration site. Chapter 6 presents the procedure used to estimate the vicarious calibration gain 
adjustment factors for bands 1-6 using the water-leaving radiances from the Marine Optical Buoy (MOBY) 
offshore of Lanai. Chapter 7 provides the adjustments to the coccolithophore flag algorithm which were required 
for improved performance over the prelaunch version. Chapter 8 is an overview of the numerous modifications to 
the atmospheric correction algorithm that have been implemented. Chapter 9 describes the methodology used 
to remove artifacts of sun glint contamination for portions of the imagery outside the sun glint mask. Finally, 
Chapter 10 explains a modification to the ozone interpolation method to account for actual time differences 
between the SeaWiFS and Total Ozone Mapping Spectrometer (TOMS) orbits. 


PROLOGUE 

The Sea- viewing Wide Field-of-view Sensor (SeaWiFS) 
Project Calibration and Validation Team (CVT) is re- 
sponsible for the overall quality of the data products and 
for verifying the processing code. The prelaunch quality 
control strategy was outlined in Volume 38 of the Sea- 
WiFS Technical Report Series (Prelaunch). Since Sea- 
WiFS began routine data processing in September 1997, 
the CVT has constantly worked to resolve data quality 
issues and improve on the initial data evaluation method- 
ologies. These evaluations have resulted in three major 
reprocessings of the entire data set (February 1998, Au- 
gust 1998, and May 2000). Each reprocessing addressed 
the data quality issues that could be identified up to the 
time of each reprocessing. 

The number of chapters (21) needed to document this 
extensive work in the SeaWiFS Postlaunch Technical Re- 
port Series requires three volumes: Volumes 9, 10, and 11. 
The chapters describe the various data quality issues, anal- 
yses, and algorithm improvements that have been devel- 
oped through the third reprocessing. The data evaluations 
after the third reprocessing indicate that the data prod- 
ucts, in most situations including Case-2 waters, are within 
the prelaunch Case-1 water accuracy goals. Nonetheless, it 
is expected that other improvements and new geophysical 
data products will be developed in the future which will 
require additional reprocessings. The SeaWiFS Project 


Office will remain dedicated to providing better products 
and to the documentation of future analysis and algorithm 
improvement studies. 

A short synopsis of each chapter in this volume is given 
below. 

j. SeaWiFS Postlaunch 

Calibration and Validation Overview 

Since launch in August 1997, the SeaWiFS Project Of- 
fice has worked diligently to improve all aspects of the 
mission, including refinements in mission operations, nav- 
igation accuracy, data processing efficiency, user and field 
operations support, and data product quality. Data prod- 
uct quality is the responsibility of the CVT. The work of 
the CVT prior to launch is largely documented in the Sea- 
WiFS Technical Report Series (Prelaunch) and other pub- 
lished documents, including journal articles. Once opera- 
tional SeaWiFS data processing began in September 1997, 
the CVT has made many improvements in the prelaunch 
processing algorithms as a result of innumerable analyses 
and evaluations. This three- volume set in the SeaWiFS 
Postlaunch Technical Report Series documents the major 
improvements and analyses that have been completed up 
to the third reprocessing in May 2000. These improve- 
ments have been incremental and have previously spawned 
reprocessings in January and August 1998. This chapter 
provides a review of the CVT activities and procedures. 
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2 . Cloud-Top Radiance Analysis for SeaWiFS 

Bilinear Gain Knee Calibration 

The distribution of cloud-top radiances in the vicinity 
of the knees of the SeaWiFS bilinear gains are used to de- 
termine the saturation counts for the ocean detectors more 
precisely than is possible based on prelaunch calibration 
data alone. Although the adjusted knee counts have little 
effect on the ocean data, they should significantly improve 
the determination of the radiances of thin dust and clouds 
in the vicinity of, and above, the knees. 

3. Lunar Data Analysis for SeaWiFS Calibration 

The SeaWiFS CVT uses monthly lunar calibrations to 
monitor the long-term stability of the radiometric calibra- 
tion of SeaWiFS. The time series of lunar observations are 
used to compute time correction factors for the SeaWiFS 
bands, which are incorporated into the calibration table. 
The response of band 7 has decreased by 3.4%, the re- 
sponse of band 8 has decreased by 9.9%, and the response 
of bands 1-6 have decreased by less than 1%, over the 
course of the mission. The calibration table will be up- 
dated as required by the on-orbit performance of the in- 
strument. 

4. Solar Data Analysis for SeaWiFS Calibration 

Daily solar calibrations are used to monitor the short- 
term variability in the radiometric response of SeaWiFS in 
between the monthly lunar calibrations. The time series 
of solar observations show the effects of both the degrada- 
tion of the reflectance of the solar diffuser and the change 
in response of the instrument. The ratio of the band 7 
time series to the band 8 time series, which minimizes the 
effects of the diffuser degradation, shows that the change 
in response of bands 7 and 8 in the solar data tracks the 
change observed in the lunar data. 

5. Vicarious Calibration of SeaWiFS Band 7 

A technique for vicariously calibrating the gain in the 
SeaWiFS band 7 is described in this chapter, as well as 
the details of its implementation. The technique has been 
found to provide a consistent estimate of the gain over the 
life of the SeaWiFS Project. 

6. MOBY Data Analysis for the Vicarious 

Calibration of SeaWiFS Bands 1-6 

The CVT performs the vicarious calibration of Sea- 
WiFS by comparing normalized water-leaving radiances 
retrieved from SeaWiFS imagery with contemporaneous 
measurements of normalized water-leaving radiances from 
the Marine Optical Buoy (MOBY). This is a system-level 
calibration that incorporates the performance of both the 
sensor and the atmospheric correction algorithm. The vi- 
carious gains for bands 1-6 depart from unity over a range 
of approximately 1-4%, depending on the band, which is 
consistent with the uncertainty in the prelaunch calibra- 
tion of SeaWiFS. 


7. The Updated SeaWiFS 
Coccolithophore Algorithm 

A coccolithophore identification algorithm has been 
adapted to work reasonably well for SeaWiFS data con- 
taining coccolithophore blooms at a number of sites in the 
Bering Sea and the North Atlantic. Changes were made 
in the initial algorithm thresholds that resulted in consid- 
erable improvements in the detection of coccolithophores 
for some cases. Future studies should include additional 
conditions which were considered in the formulation of 
the original coccolithophore algorithm, but not considered 
here— red tides, sediments, whitings, and haze. Finally, a 
future study may want to re-examine the predictors that 
are employed to detect the coccolithophore signal. The 
SeaWiFS 443, 510, and 555 nm bands are slightly different 
than the CZCS bands used to derive the current predic- 
tors. In addition, SeaWiFS has bands at 412 and 490 nm, 
which CZCS did not have, and these might be useful in de- 
riving a better set of predictors to detect coccolithophores 
and other water types. 

8. The SeaWiFS Atmospheric 
Correction Algorithm Updates 

Modifications to the atmospheric correction algorithm 
for the third SeaWiFS reprocessing are described in this 
chapter. The updates include changes to the aerosol look- 
up tables, the atmospheric diffuse transmittance tables, 
the ocean whitecap computations, and the implementation 
of new Rayleigh radiance tables generated with a variety 
of ocean surface wind speeds. In addition, computation 
of a new SeaWiFS atmospheric product, the Angstrom 
exponent, is described. These modifications significantly 
improve the SeaWiFS retrieval results. 

9. Correction of the Sun Glint 
Contamination on the SeaWiFS 

Aerosol Optical Thickness Retrievals 

For ocean color remote sensing, the measurement of ra- 
diances affected by sun glint has to be avoided or masked 
out. SeaWiFS has a capability of operationally tilting the 
sensor 20° away from nadir to minimize sun glint contami- 
nation, however, sun glint is still a factor near the subsolar 
point. In this chapter, results are presented which quan- 
tify the effect of sun glint contamination on the retrievals 
of atmospheric and bio-optical oceanic products. It was 
found that, although the sun glint contamination has a 
minor effect on the retrieved bio-optical oceanic products, 
the effect on the retrieved atmospheric products (e.g., aero- 
sol optical thickness) is significant. A sun glint correction 
scheme is described, which was implemented in the Sea- 
WiFS data processing. It was found that the sun glint 
correction significantly improves the derived atmospheric 
products in the vicinity of the subsolar point. 
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10. Modifications to the TOMS 

Ozone Ancillary Data Interpolation 

The method for interpolating the Total Ozone Mapping 
Spectrometer (TOMS) ozone data used for level -2 process- 
ing was changed as a result of new information about the 
generation of the files by the TOMS Project. The new 
method and comparison of the old and new results are de- 
scribed in this report. 
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Chapter 1 


SeaWiFS Postlaunch Calibration and 
Validation Overview 

Charles R. McClain 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 

Abstract 

Since launch in August 1997, the SeaWiFS Project Office has worked diligently to improve all aspects of the 
mission, including refinements in mission operations, navigation accuracy, data processing efficiency, user and 
field operations support, and data product quality. Data product quality is the responsibility of the CVT. The 
work of the CVT prior to launch is largely documented in the SeaWiFS Technical Report Senes (Prelaunch) 
and other published documents, including journal articles. Once operational SeaWiFS data processing began 
in September 1997, the CVT has made many improvements in the prelaunch processing algorithms as a result 
of innumerable analyses and evaluations. This three- volume set in the SeaWiFS Postlaunch Technical Report 
Series documents the major improvements and analyses that have been completed up to the third reprocessing 
in May 2000. These improvements have been incremental and have previously spawned reprocessings in January 
and August 1998. This chapter provides a review of the CVT activities and procedures. 


1.1 INTRODUCTION 

The SeaWiFS calibration and validation program has 
four basic areas of responsibility: 

1) Providing the processing algorithms; 

2) Verifying and documenting the level conversion pro- 
cessing code, i.e., level-0 through level-3; 

3) Tracking sensor performance and calibration from 
fabrication throughout the lifetime of the mission; 
and 

4) Verifying the quality of the data products. 

These responsibilities encompass a broad range of scien- 
tific and technological areas of expertise. Given the orig- 
inal 28-month launch schedule, the SeaWiFS Project was 
hard pressed to define and implement a viable capability 
for mission operations, data capture, calibration and vali- 
dation, and data processing. 

For the CVT, the principal group that was originally 
identified by NASA management to assist the SeaWiFS 
Project was the Moderate Resolution Imaging Spectrora- 
diometer (MODIS) Oceans Team. The SeaWiFS Project 
directly funded members of that team (H. Gordon, R. 
Evans, D. Clark, and K. Carder) to accelerate their re- 
spective activities and to assist the SeaWiFS Project in 
preparing for launch. 

The CVT also sought the participation and assistance 
of the ocean color community, fully recognizing that suc- 


cessful fulfillment of its objectives required the participa- 
tion of a broader community. In many cases, this assis- 
tance was funded under grants and contracts, but most 
often, it was provided by individual researchers on a vol- 
untary basis. As part of this process, the CVT identified 
areas where the community infrastructure and capability 
was lacking or nonexistent for supporting a global ocean 
color mission. 

To address the identified deficiencies, the CVT and the 
SeaWiFS Project Office initiated or accelerated a number 
of efforts: 

1) Measurement protocols; 

2) Instrument intercalibrations organized as SeaWiFS 
Intercalibration Round- Robin Experiments, (SIR- 
REXs); 

3) Bio-optical data archiving within the SeaWiFS Bio- 
optical Archive Storage System, (SeaBASS); 

4) Prelaunch Bio-optical Algorithm Working Group! ; 

5) Advanced laboratory and m situ measurement tech- 
nology; 


f After launch, the Sensor Intercomparison and Merger for 
Biological and Interdisciplinary Oceanic Studies (SIMBIOS) 
program assumed many of the responsibilities related to al- 
gorithm development. 
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6) SeaWiFS Technical Report Series f and; 

7) Imagery visualization and analysis with the Sea- 
WiFS Data Analysis System (SeaDAS)J 

Over the course of the extended prelaunch- and operational 
phases of the project, the CVT conducted and participated 
in many activities which are summarized in Table 1 and are 
reviewed in Hooker and McClain (2000). 

The revised set (the third reprocessing) of archive geo- 
physical data parameters and the quality control masks 
and flags are listed in Table 2. Level-la data are in sensor 
counts and is navigated, but requires a separate transfor- 
mation provided by the CVT to convert counts to radi- 
ance. The level *2 and -3 products are in geophysical units, 
e.g., milligrams per cubic meter (mgm“ 3 ). It is expected 
that the archive product list will grow as the ocean com- 
munity defines new products and algorithms and as other 
Earth science disciplines begin using SeaWiFS data for at- 
mospheric and terrestrial applications. Indeed, the Sea- 
WiFS Project is actively encouraging the development of 
new products and applications and has worked with other 
groups to generate non-archive evaluation products, such 
as the normalized difference vegetation index (NDVI) over 
land. The SeaWiFS Project, however, does not have the 
expertise nor the resources to verify non-oceanic products. 

1.2 DATA PROCESSING ALGORITHMS 

Data processing algorithms include those for sensor cal- 
ibration, stray light, atmospheric correction, bio-optical 
properties, masks, flags, and level -3 binning. There is a 
substantial amount of literature available on these topics 
for which Table 1 provides many of the references. It is 
not sufficient to simply verify and document individual 
algorithms — the organization, links, and sequences of data 
manipulations within the data processing system must also 
be described (Darzi 1998). A major component of all this 
is the creation of the file specifications and input-output 
routines. Much of this was undertaken by the CVT. In 
fact, SeaWiFS was one of the first NASA projects to adopt 
the hierarchical data format (HDF) and worked with the 
National Center for Supercomputing Applications (NCSA) 
during the HDF development. 

The prelaunch verification of the processing was greatly 
facilitated by creating simulated data (Gregg et al. 1993 
and 1994) by the Mission Operations Element and routine 


f The SeaWiFS Technical Report Series (pre- and postlaunch) 
are funded by the SeaWiFS Project Office. The majority 
of volumes are authored by members of the CVT and its 
collaborators. 

X The SeaDAS interactive image processing software currently 
supports processing from the Coastal Zone Color Scanner 
(CZCS), SeaWiFS, the Ocean Color and Temperature Scan- 
ner (OCTS), and the Marine Optical Spectroradiometer 
(MOS). In addition, it supports the display of MODIS ocean 
products. 


end-to-end processing coordinated by the Data Process- 
ing Element. The CVT worked closely with the algorithm 
providers, especially the University of Miami collaborators 
(H. Gordon, R. Evans, J. Brown, and others), who pro- 
vided much of the level -2 and -3 processing code. As a re- 
sult of the extensive prelaunch testing, data flowed through 
the processing system and to the GSFC Distributed Active 
Archive Center (DAAC) within hours of the first OrbView- 
II transmission to the receiving station at the NASA Wal- 
lops Flight Facility. 

During the 120 day postlaunch data acceptance period, 
the CVT performed extensive analyses to verify that the 
acceptance criteria (Hooker and McClain 2000) were sat- 
isfied and presented its findings at a miniworkshop with 
members of the MODIS Oceans Team. Corrections to 
the operational algorithms during the first four months 
culminated in the first reprocessing in January-February 
1998. More substantial changes, including incorporating 
time-dependent adjustments to the sensor calibration, ne- 
cessitated the second (August-September 1998) and third 
(May 2000) reprocessings. Robinson et al. (2000) summa- 
rized the primary changes for each reprocessing. Prepa- 
rations for the third reprocessing included two miniwork- 
shops where algorithm issues were reviewed and evaluation 
strategies were defined. The results of analyses conducted 
in preparation for the second and third reprocessings were 
posted on project web sites for consideration by the user 
community. Finally, for the third reprocessing, all level -1, 
-2, and -3 processing codes were rewritten and stream- 
lined. 

1.3 SENSOR CALIBRATION 

Table 1 is a list of the various activities and accomplish- 
ments of the CVT including documents dealing with the 
SeaWiFS calibration and characterization and processing 
algorithms. In the postlaunch phase, the CVT has pursued 
a variety of analyses for tracking the on-orbit behavior of 
the instrument, which are included in the following list. 

A. Prelaunch 

1. Laboratory sensor characterization and calibra- 
tion 

2. Solar calibration for a transfer- to-orbit compar- 
ison 

B. Postlaunch 

1. Operational Adjustments 

a. Lunar calibration (monthly) for time depen- 
dence correction 

b. Solar calibration (daily) for time dependence 
(bands 7 and 8) and fine resolution check of 
lunar correction 

c. Open ocean e(765, 865) analysis for vicarious 
calibration (band 7) 

d. MOBY Lwn time series for vicarious cali- 
bration (bands 1-6) 
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Table 1. Specific activities and accomplishments of the SeaWiFS CVT and collaborators. The term “Vol.” indicates 
a SeaWiFS Technical Report Series Prelaunch volume (Appendix A), and the term “Vol. ” indicates a Posfciaunch 
volume (Appendix A). The CVT web site is located at http://seawifs.gsfc.nasa.gov/~grey/calval.html. 


Activity or Accomplishment 

Reference 

Round-Robins and Measurement Protocols 

Seven instrument calibration round-robins 

Vols. 14, 16, 34, 37, and VOLS. 7 

Joint US-Japan SeaWiFS-OCTS 

Johnson et al., 1997 

prelaunch cross-calibration 

Johnson et al., 1997 

Data collection protocols 

Vols. 5 and 25 (revision) 

Data analysis round-robin 

Vol. 26 

High performance liquid chromatography (HPLC) round-robin 

Vol. 14 

Tower shading effects experiment 

Zibordi et al. 1999 

Two measurement protocol experiments 

Vols. 3 and 8 

Advanced Instrumentation and Community Support 

SeaWiFS SWG in situ radiometer calibration support 

CHORS, 1992-1995 

SeaWiFS Transfer Radiometer (SXR) 

Vol. 1 

SeaWiFS Quality Monitor (SQM) 

Johnson et al. 1998; Hooker and Aiken 1998 

SeaPRISM 1 development with the JRC 2 

Vol. 13 

NASA Environmental Sampling System development 

Hooker and Maritorena 2000 

| SeaWiFS Sensor Calibration and Characterization } 

Two prelaunch solar-based calibration experiments 

Vols. 19 and 27 

Prelaunch acceptance report 

Vol. 22 

Prelaunch calibration reports 

Vol. 23 and Vol. 4 

Stray light description 

Vol. 31 

Solar diffuser design 

Vol. 39 

Calibration temperature dependence 

Vol. 40 

Postlaunch data acceptance evaluation 

December 1997 

Science quality data certification 

McClain et al. 1998 

Lunar and solar data analysis procedures 

Barnes et al. 1999a 

Prelaunch solar calibration transfer to orbit analysis 

Vol. 5 

1 Co-Funded MODIS Oceans Team Activities 

Atmospheric corrections 

Gordon 

Semi-analytical chlorophyll algorithm 

Carder 

In situ vicarious calibration 

Clark 

Three operational MOBYs and one shipboard spectrometer 

Clark et al. 1997 

Support facility in Honolulu and routine deployments 

July 1997-present 

Two Marine Optical Characterization Experiments (MOCE) 

Clark 

Postlaunch data acceptance workshop 


SeaWiFS postlaunch initialization cruise 

Clark and Gordon 

Additional Supported Investigations 

Atmospheric correction studies 

Vols. 13, 19, and 27 

Bermuda Atlantic Time Series (BATS) observations 

Siegel 

California Cooperative Fisheries (CalCOFI) cruises 

Mitchell 

Plymouth Marine Bio-Optical Data buoy (PlyMBODy) 

Vol. 33 

Nine Atlantic Meridional Transect (AMT) cruises with PML 3 

Vol. 35, VOL. 2, and Aiken et al. 2000 

Portable optical measurements laboratory 


CTD and Niskin bottle rosette 


Bio-optical algorithm development and evaluation 

O’Reilly 

Acqua Alta Oceanographic Tower (AAOT) observations 

Zibordi 

Productivity des Systemes Oceaniques Pelagiques (PROSOPE) 4 

Claustre and Morel 
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Table 1. (cont.) Specific activities and accomplishments of the SeaWiFS CVT and collaborators. The term “VOL.” 
indicates a SeaWiFS Technical Report Series Prelaunch volume (Appendix A), and the term “VOL. ” indicates a Post- 
launch volume (Appendix A). The CVT web site is located at http : //seawif s . gsf c . nasa . gov/~grey/ calval . html. 


Activity or Accomplishment 

Reference 

Algorithm Development 

Workshops (Meeting summaries published) 

Vols. 18, 24, 36, and 43 

Seven prelaunch bio-optical algorithm and 
protocols workshops 

Absorption measurement workshop at SIO 5 

Mitchell 

One calibration workshop 

Two atmospheric correction workshops 

(One was a joint meeting with MODIS Project) 
SeaWiFS Bio-optical Algorithm Mini-workshop (SeaBAM) 
SeaWiFS postlaunch data acceptance workshop 
Two SeaWiFS algorithm evaluation mini-workshops 

Atmospheric correction algorithm 

Gordon and Wang 1994a, and Gordon 1995 

Bio-optical data archive SeaBASS 

Vol. 20 

Data quality masks and flag algorithms 

Vol. 28 

Out-of-band effects and correction scheme 

VOLS. 28, 39, 40, and 41 

Initial level -3 binning algorithm 

Vol. 32 

K( 490) algorithm 

VOL. 41 and VOL. 11 

SeaWiFS stray-light correction algorithm 

Vol. 41 

CZCS pigment and chlorophyll a algorithms 

O’Reilly et al. 1998 and VOL. 11 

CVT image gallery 

http : //calval-2 . gsf c . nasa . gov/calval/ 

Quality Control (QC) Software 

Operational QC Software 

Vol. 38 

Level-1, -2, and -3 quality assurance 
Ancillary data quality assurance (winds, pressure, ozone) 
Derived product evaluation 
Sensor engineering telemetry data tracking 
(Presently handled by Mission Operations) 

Calibration evaluation (lunar, solar, vicarious) 

Operational Processing Software 

Format specifications and code for all products 

http : //seawif s . gsf c . nasa . gov/SEAWIFS/ 
SOFTWARE/SOFTWARE . html 
#Product Specifications 

Documentation on final processing flow completed 

Darzi 1998 

Level-1, -2 and -3 code verification 

(At-launch and first through third reprocessings) 

SeaWiFS Data Analysis System 6 (SeaDAS) 

Three training classes 

Summer 1994 

Preliminary version 

Summer 1994 

Interactive Data Language (IDL) licenses (45 total) 

Distributed to the Science Working Group 

One training class 

May 1995 

Version 1 

Summer 1995 

Sun workstation delivered to Shirshov Institute (Moscow) 

1995 

One training class 

April 1996 

Version 2 

May 1996 

Version 3 

September 1997 

Seven training classes 

November 1997 
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Table 1 . (cont.) Specific activities and accomplishments of the SeaWiFS CVT and collaborators. The term “Vol.” 
indicates a SeaWiFS Technical Report Series Prelaunch volume (Appendix A), and the term “VOL. ” indicates a Post- 
launch volume (Appendix A). The CVT web site is located at http://seawifs.gsfc.nasa.gov/~grey/calval.html. 


Activity or Accomplishment 

Reference 


SeaDAS (cont.) 

Version 3.1 

February 1998 

Two training classes 

July 1998 

Version 3.2 

October 1998 

Version 3.3 

April 1999 

Linux version 

May 1999 

Version 4 

May 2000 


1. SeaPRISM is the SeaWiFS Photometer Revision for Incident Surface Measurement. 

2. JRC is the Joint Research Centre (Ispra, Italy). 

3. PML is the Plymouth Marine Laboratory (Plymouth, England). 

4. The French meaning for PROSOPE translates to productivity of pelagic oceanic systems. 

5. SIO is the Scripps Institution of Oceanography. 

6. SeaDAS is primarily funded by the NASA Oceanography Program, however, the SeaWiFS Project is presently providing 
the funding needed for equipment upgrades and system administration assistance. The SeaDAS staff rely largely on the 
CVT for processing code. 


Table 2. SeaWiFS archived atmospheric and ocean surface products for the third reprocessing. The QC masks 
and flags are listed as well because some are used as exclusion criteria for the level -3 binning. The Tilt and Sensor 
Engineering Limits flags are applied line-by-line and are not represented in the level -2 products as graphic overlays. 


Level 

Product Type 

Product Name 

Level -1 

Geolocated Sensor Countsf 


Level -2 

Ocean Products 

Lwn( 412, 443, 490, 510, 555, and 670 nm) accuracy goal: ±5% 
Chlorophyll a accuracy goal: ±35% in Case-1 water 
K( 490) 


Atmospheric Products 

7(510, 865) (Angstrom exponent) 
Ta (865) 

«(765, 865) 


QC Masks 

LandJ 

Cloud and icej 
Sun glintt 

Atmospheric correction failure or invalid data} 
High L t i . 


QC Flags 

Large solar zenith angle (increased to 75°) J 
Large satellite zenith anglej 
Negative Lw (bands 1-5) § 

Stray light J 
CoccolithophoreJ 
Low Lwn(5 55)f 

Outside chlorophyll a algorithm range (0-64.0) 

Missing ancillary data 

Turbid Case-2 water 

Shallow water 

High r a (865) 

Tilt underway:): 

Absorbing aerosol^ 

Trichodesmium 
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Table 2. (cont.) SeaWiFS archived atmospheric and ocean surface products for the third reprocessing. The 
QC masks and flags are listed as well because some are used as exclusion criteria for the level -3 binning. The 
Tilt and Sensor Engineering Limits flags are applied line-by-line and are not represented in the level -2 products 
as graphic overlays. 


Level 

Product Type 

Product Name 

Level -2 

QC Flags 

Chlorophyll algorithm failure} 

Maximum number of NIR iterations exceeded} 
Sun glint correction 
L r > Li 

Atmospheric correction warning} 

Level “3 

Binned Products 

All level *2 fields 

C a 'K [ratio of chlorophyll a/K (490)] 


f CVT provides calibration tables to convert counts to radiance separately. 

} Masks and flags used as exclusion criteria in the generation of leveT3 binned products. 
§ Values set to zero for binning. 


C. Product Evaluations 

1. Global clear- water time series (8-day binned 
data) using Lwn (bands 1-6) e(765,865), 
r a (865), chlorophyll a, number of clear-water 
bins 

2. L t - L r < 0 and negative Lw analyses (statistics 
and global distributions) 

3. r 0 (865) comparisons with Aerosol Robotic Net- 
work (AERONET) data and SeaWiFS 

4. Lw and Lwn comparisons (in situ versus Sea- 
WiFS) 

5. Chlorophyll a and #(490) comparisons (in situ 
versus SeaWiFS) 

6. E s (surface irradiance) comparisons (in situ ver- 
sus theoretical clear-sky values) 

7. Global time series of the mean cloud albedo at 
865 nm (inconclusive and discontinued) 

8. Earth curvature effects on L r estimates 

9. f/Q corrections (bidirectional reflectance) to 
Lw values (Morel and Gentili 1996) 

10. Effects of alternative solar spectra, Fo(A), on vi- 
carious calibrations and derived products 

11. Band 8 calibration evaluation using clear sky, 
low chlorophyll region along Antarctic coast 

Of the postlaunch analyses, all are described in detail in 
various chapters of this volume with the exception of the 
cloud albedo time series. The cloud albedo time series 
proved to be very noisy; it was dropped after the lunar 
and solar calibrations were shown to be very robust for 
tracking sensor degradation. 

The overall scheme for sensor calibration is shown in 
Fig. 1; it illustrates how the various prelaunch and post- 
launch calibrations are connected. The lunar calibration 
is used to remove any time dependence in the sensitiv- 
ity (Eplee and Barnes 2000); the solar calibration is used 


as a cross check of the lunar measurements for bands 7 
and 8 (Eplee et al. 2000). Because there is no accurate 
method for vicariously calibrating band 8 at this time, the 
prelaunch calibration is assumed. This assumption is ver- 
ified to within ±2% by the solar calibration transfer-to- 
orbit results (Barnes et al. 1999b). These results are sup- 
ported by an analysis of clear sky (cloud albedo thresh- 
old set at 0.35), low chlorophyll waters, i.e., no NIR re- 
flectance, along the Antarctic coast near 60°S,25°E during 
November 1997 through January 1998. Under these condi- 
tions, Lt(865) always exceeded L r (86 5), i.e., band 8 is not 
under calibrated, and the average r a (865) value was around 
0.01, similar to the minimum values observed at the Mc- 
Murdo research station in the Antarctic during the same 
time, i.e., band 8 is not significantly overcalibrated. 

After the time dependencies are removed, the prelaunch 
calibration coefficients (Johnson et al. 1999) for bands 1-6 
are adjusted using MOBY matchups to minimize the aver- 
age differences between the buoy and SeaWiFS normalized 
water- leaving radiances (Eplee and McClain 2000). The 
band 7 (765 nm) calibration is adjusted so the atmospheric 
correction parameter at 765 and 865 nm, e( 765, 865), is 
near the expected value at the MOBY site (Robinson and 
Wang 2000). Checks of the results of this process include 
analyses of the eight-day global binned products to ver- 
ify that no unexpected trends are occurring. For example, 
the second global reprocessing was initiated when it was 
realized that bands 7 and 8 were degrading. This degrada- 
tion resulted in a steady increase in the e(765, 865) values 
with a commensurate gradual decrease in the global mean 
clear-water radiances. 

Match-up analyses using independent in situ data are 
also used for verification (Bailey et al. 2000). These analy- 
ses show that the SeaWiFS clear-water radiances after the 
second reprocessing compare very well with the in situ val- 
ues suggesting that the sensor calibration is correct. The 
comparisons, however, over regions of high chlorophyll and 
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Fig. 1 . Schematic of the SeaWiFS calibration chronology and methods. 


turbid water show that the SeaWiFS water-leaving radi- 
ances are low, particularly at 412 and 443 nm, which has 
been attributed to the assumption of zero water-leaving 
radiance at 765 and 865 nm in the atmospheric correction 
algorithm. One of the primary reasons for initiating the 
global third reprocessing is to address this problem. Below 
is an outline of the entire calibration verification procedure. 

A. Determine temporal degradation lunar calibration data 

1. Degradations in bands 1, 2, 5, 6, 7, and 8 observed 

B. Determine nominal Lwn values from MOBY 

1. Use the same values as used in MOBY-SeaWiFS 
match-up data set 

C. Set band 7 calibration correction factor 

1. Assume band 8 prelaunch calibration 

a. Calibration correction factor = 1.0 

b. Time-dependent degradation correction from lu- 
nar data 

2. Adjust band 7 calibration factor to get appropriate 
mean e(765, 865) value 

a. Time-dependent degradation correction from lu- 
nar data 

b. Use MOBY LAC time series 

D. Analyze MOBY Lwn data for bands 1-6 
1. Use theoretical E s values 


2. Check for trends in time series 

3. Use MOBY LAC time series 

4. Apply exclusion criteria to match-up data set 

5. Adjust calibrations to minimize SeaWiFS-MOBY 
Lwn differences for valid match-up data 

a. Compute geometric mean of match-up ratios 

b. Generate final calibration adjustment factors 

E. Check global clear- water Lwn values from eight-day 
composite time series and compare them with nominal 
MOBY clear-water values 

F. Evaluate mask and flag performance (e.g., cloud and 
coccolithophore) 

G. Run match-up analyses on other ship and buoy data 
for verification of clear-water Lwn values 

1.4 PRODUCT VALIDATION 

Product validation consists of match-up analyses and 
real-time quality control. McClain et al. (1996) outline 
the various analyses that were envisioned prior to launch. 
Since launch, these analyses have been substantially ex- 
panded and refined as is outlined in this volume. All global 
area coverage (GAC) and local area coverage (LAC), in- 
cluding high resolution picture transmission (HRPT) sta- 
tion data products, are passed through automated checks 
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and visual inspections before being approved for archiv- 
ing at the GSFC DAAC. Not every file can be reviewed 
visually, but those that trigger a flag from the automated 
checks are inspected. 

Another aid in diagnosing problems is the routine gen- 
eration of quality control products which are not archived 
and are automatically purged after a certain period of time. 
The quality control products up to the third reprocessing 
were: 

a) Zonal wind; 

b) Ozone; 

c) Sensor azimuth angle; 

d) L r ( 443); 

e) Second aerosol model; 

f) Merdional wind; 

g) Solar zenith angle; 

h) Lvw(670); 

i) L s (865); 

j) Surface pressure; 

k) Albedo at 865 nm; 

l) M765); 

m) L/(865); 

n) Precipitable water; 

o) Sensor zenith angle; 

p) M443); 

q) First aerosol model; and 

r) e(765, 865) flag, i.e., pixels having values outside the 
range of valid e(765,865) values. 

These products were discontinued once the third reprocess- 
ing commenced, as their routine generation was no longer 
necessary for product evaluation. 

The other aspect of quality control is the definition of 
quality control masks and flags. These are listed in Ta- 
ble 2 along with the level-3 exclusion criteria. Pixels that 
are masked are not processed to level -2. Pixels that are 
flagged are processed, but are not necessarily included in 
the level -3 products. Beginning with the third reprocess- 
ing, suites of quality control products and masks and flags 
will be refined at each reprocessing. It should be noted 
that the present suit of masks and flags are designed to 
optimize the accuracy of the chlorophyll a product. In the 
future, as suites of products for atmospheric, terrestrial, 
and other ocean parameters are defined, each will need its 
own set of masks and flags. 

Finally, the CVT works closely with the SIMBIOS 
Project (McClain and Fargion 1999), which is co-located 
with the SeaWiFS Project, to provide cruise support so 
as to optimize the collection of in situ data with respect 
to the SeaWiFS coverage. This support includes overpass 
predictions in advance to assist researchers planning ship 
tracks and station times. The project provides real-time 


data products, which can be tailored to the researchers * 
needs via the SIMBIOS web sitef. These products can be 
electronically mailed to the ship or to a point of contact. 
In addition, if the ship locations are known a week or more 
in advance of the satellite overpass, onboard LAC data can 
be scheduled to ensure that high resolution data over the 
ship will be available for match-up analyses. 

During the first two years of operations, the project 
supported 125 field studies. In addition, LAC data are 
routinely scheduled over time series sites, such as MOBY, 
the AAOT site (offshore of Venice, Italy); the BATS site; 
the Hawaii Ocean Time-series (HOT) site; certain Trop- 
ical Ocean Global Atmosphere (TOGA) program Tropi- 
cal Atmosphere-Ocean (TAO) moorings; and others. LAC 
recorder space that is not used to cover validation targets is 
used over default regions of interest, such as the Galapagos 
Islands. Figure 2 provides a typical LAC data collection 
summary. The SIMBIOS program supports a number of in 
situ data collection activities including cruise and mooring 
bio-optical data and aerosol optical thickness from a num- 
ber of AERONET sites. The match-up methodology for 
atmospheric parameters is outlined in Wang et al. (2000). 

1.5 SUMMARY 

Since the inception of the SeaWiFS Project in the early 
1990s, the CVT has labored diligently to ensure that a 
comprehensive and effective calibration and validation pro- 
gram was in place by launch. The four-year launch delay 
allowed time for many capabilities to be realized in time 
for launch. The overall philosophy of the CVT has been 
to involve the research community as partners and to ini- 
tiate activities which develop community infrastructure. 
In order to gain a better appreciation for the challenges 
of field data collection, the CVT has developed its own 
field program dedicated to improving measurement accu- 
racy and collecting high quality bio-optical data. Dr. Stan- 
ford Hooker leads this effort and has been actively involved 
in the AMT program (Robins et al. 1996 and Aiken et al. 

1998) , the PROSOPE cruise, and a variety of measurement 
protocol experiments on the AAOT platform (Hooker et al. 

1999) . 

With each reprocessing, significant improvements in 
the data products have been achieved. After the second 
reprocessing, problems with low or negative water-leaving 
radiances persisted in certain situations. As a result, im- 
provements have been made in a number of algorithms in 
preparation for the third reprocessing: 

1. Sensor degradation correction (Eplee and Barnes 

2000 ); 

2. Bilinear gain knee offset adjustments (Eplee and 
Patt 2000); 

3. Improved ozone interpolation scheme (Ainsworth 
and Patt 2000); 


f SIMBIOS URL: http://simbios.gsfc.nasa.gov 
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Fig. 2. A typical daily SeaWiFS onboard LAC data collection schedule. 


4. Sun glint correction algorithm (Wang and Bailey 

2000 ); 

5. Surface whitecap correction (Robinson et al. 2000); 

6. Various atmospheric correction algorithm improve- 
ments (Wang 2000 and Siegel et al. 2000); 

7. Absorbing aerosol detection (Hsu et al. 2000); 

8. Out-of-band corrections to L w values (Wang et al. 
2000 ); 

9. A (490) algorithm (Mueller 2000); and 

10. Improved chlorophyll a algorithms (O’Reilly et al. 

2000). 

These improvements were incorporated in the third repro- 
cessing. 

Other potential improvements were evaluated, but not 


incorporated. For instance, the Rayleigh radiance model 
was compared to other models, some with Earth curva- 
ture effects included. No significant differences were found 
for solar zenith angles that were less than 70°. In addi- 
tion, the bidirectional reflectance algorithm of Morel and 
Gentili (1996) was tested, but the use of the algorithm 
was deferred until certain improvements could be imple- 
mented. In addition, a different solar spectrum was tried, 
but it yielded almost identical results as the current spec- 
trum, because the vicarious calibration compensates for 
the differences. 

As further improvements and new products are defined, 
annual reprocessings are anticipated; the CVT will con- 
tinue to work with the Data Processing Element on how 
processing efficiency can be maximized as the SeaWiFS 
data set grows. 
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Chapter 2 


Cloud- Top Radiance Analysis for SeaWiFS 
Bilinear Gain Knee Calibration 

Robert E. Eplee, Jr., and Frederick S. Patt 
SAIC General Sciences Corporation 
Beltsville, Maryland 

Abstract 

The distribution of cloud-top radiances in the vicinity of the knees of the SeaWiFS bilinear gains are used to 
determine the saturation counts for the ocean detectors more precisely than is possible based on prelaunch 
calibration data alone. Although the adjusted knee counts have little effect on the ocean data, they should 
significantly improve the determination of the radiances of thin dust and clouds in the vicinity of, and above, 
the knees. 


2.1 INTRODUCTION 

An examination of histograms of cloud-top radiances 
observed by SeaWiFS shows anomalous distributions of ra- 
diances in the vicinity of the knee of the bilinear response 
for each SeaWiFS band. These anomalous distributions 
are observed in HRPT, LAC, and GAC images. Typical 
histograms for HRPT images, which are plotted for each 
band in Fig. 3, show unexpected discontinuities in the dis- 
tribution of radiances in the vicinity of the knees: there 
are fewer pixels at the knee radiances than expected from 
the number of points with radiances below and above the 
knees. These discontinuities in the radiance distribution 
would give rise to errors in the retrieved radiances of thin 
dust and thin clouds, as measured by SeaWiFS. 

Each SeaWiFS band has three ocean detectors with 
high radiometric sensitivity and low saturation radiances, 
and one cloud detector with low radiometric sensitivity and 
high saturation radiances. During the collection of ocean 
data, SeaWiFS averages the output from the four detec- 
tors within each band. The knees in the response occur 
where the ocean detectors saturate. Figure 4 shows bilin- 
ear response for bands 1-8 in the vicinity of the knees. A 
plot of the biliner response of band 1 over the full dynamic 
range of the instrument is shown in Fig. 8 of Johnson et 
al. (1999). 

The discontinuities in radiance in the vicinity of the 
knees of the bilinear response are caused by uncertainties 
in determining the output counts at which saturation oc- 
curs for the ocean detectors. The prelaunch radiometric 
calibration of SeaWiFS was performed in the laboratory 
at four radiance levels, which were either below or above 


the knee radiances (Johnson et al. 1999). Sunlight reflected 
from cloud tops provides a continuum of radiances in the 
vicinity of the knees. The CVT has undertaken an analy- 
sis of the distribution of cloud-top radiances measured by 
SeaWiFS to make a more precise determination of the satu- 
ration counts for the ocean detectors than is possible based 
on the prelaunch calibration data alone. 

2.2 DETERMINING KNEE COUNTS 

The CVT can adjust the counts at which saturation 
occurs for each detector in the SeaWiFS calibration model 
by changing the offsets used in computing the counts-to- 
radiance conversion factors from the detector gains (Darzi 
et al. 1995). The CVT performed the adjustment for a 
given band with a trial-and-error approach, where the off- 
sets for each detector were changed individually and the 
resulting radiance distribution was examined. The goal of 
these adjustments was to improve the distribution of ra- 
diances in the vicinity of the knees for each band. The 
final offsets were those which yielded the smoothest radi- 
ance distrubutions at the knees. Because SeaWiFS only 
collects ocean data for gain 1, this knee calibration can 
only be performed at this gain. 

Figure 3 also shows the radiance distributions result- 
ing from the modified knees for bands 2-8. For each of 
these bands, the distribution for the modified knees is an 
improvement over the original distribution. A correction 
could not be determined for band 1 because a one count 
adjustment of the saturated counts for this band resulted 
in a poorer radiance distribution. Figure 4 shows the lo- 
cations of the modified knees for each band. 
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Fig. 3. Histograms of radiance distributions in the vicinity of the knees in the bilinear radiometric 
response. 
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Fig. 4. SeaWiFS bilinear radiometric response in the vicinity of the knees. 
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Table 3. Knees of the original SeaWiFS bilinear gains for gain 1. 


Band 

Knee 1 

Radiance Counts 

Knee 2 

Radiance Counts 

Knee 3 

Radiance Counts 

Saturation 
Radiance Counts 


10.98 

793.23 

11.00 


11.14 

797.98 

60.69 

1002.25 


10.60 

789.48 

10.62 

1 

10.68 

792.24 

68.85 

1004.50 


8.344 

799.96 

8.356 


8.395 

782.03 

69.46 

1002.25 


7.171 

778.31 

7.174 

778.50 

7.189 

779.09 

67.10 

1002.75 

5 

5.871 


5.887 

772.48 

5.920 

773.96 

67.19 

1001.25 

6 

3.324 

762.31 

3.326 

762.70 

3.342 

763.93 

57.05 

999.75 

7 

2.361 

758.99 

2.379 

762.85 

2.389 

763.93 

43.64 

1000.25 

8 

1.695 

762.23 

1.699 

763.46 

1.704 

764.23 

34.94 

1002.75 


Table 4. Knees of the modified SeaWiFS bilinear gains for gain 1. 


Band 

Knee 1 

Radiance Counts 

Knee 2 

Radiance Counts 

Knee 3 

Radiance Counts 

Saturation 
Radiance Counts 

1 

10.98 

793.23 

11.00 

794.10 

11.14 

797.98 

60.69 

1002.25 

2 

10.61 

790.27 

10.63 

791.16 

10.71 

793.31 

68.85 

1005.50 

3 

8.344 

799.96 

8.356 

780.71 

8.404 

782.31 

69.46 

1002.50 

4 

7.171 

778.31 

7.174 

778.50 

7.203 

779.64 

67.10 

1003.25 

5 

5.877 

771.78 

5.887 

772.73 

5.920 

774.21 

67.19 

1001.50 

6 

3.330 

763.83 

3.333 

764.22 

3.352 

765.72 

57.05 

1001.50 

7 

2.363 

759.75 

2.379 

763.10 

2.391 

764.45 

43.64 

1000.75 

8 

1.697 

762.99 

1.701 

764.22 

1.709 

765.52 

34.94 

1004.00 


The change in location of the knee radiances and counts 
for each band is also shown by comparing Tables 3 and 4, 
which list the radiance and counts at the three knees before 
and after the adjustments of the offsets. The changes in 
the saturated counts at the knees are from 0.0-1. 8 counts, 
with a mean of 0.6 counts. The typical change in the ra- 
diance for a given instrument output below the knee is 
0.1%, which is within the uncertainty of the instrument 
calibration (Johnson et al. 1999). The typical change in 
the radiance above the knee is 0.8%, which again is within 
the uncertainty of the instrument calibration (Johnson et 
al. 1999). The changes in the above-the-knee radiances 
should improve the determination of dust and absorbing 
aerosol masks for SeaWiFS (Hsu 2000). The offsets for the 


modified knee locations have been incorporated into the 
SeaWiFS calibration table which will be used for the third 
reprocessing of SeaWiFS data. 

2.3 DISCUSSION 

By using the distribution of cloud-top radiances in the 
vicinity of the knees of the SeaWiFS bilinear gains, the 
CVT made more precise determinations of the knee counts 
than were possible from the prelaunch calibration data. 
While the adjusted knee counts have little effect on the 
ocean data, they should significantly improve the deter- 
mination of the radiances of thin dust and clouds in the 
vicinity of, and above, the knees. 
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Chapter 3 


Lunar Data Analysis for SeaWiFS Calibration 

Robert E. Eplee, Jr., and Robert A. Barnes 
SAIC General Sciences Corporation 
Beltsville, Maryland 

Abstract 

The SeaWiFS CVT uses monthly lunar calibrations to monitor the long-term stability of the radiometric cali- 
bration of SeaWiFS. The time series of lunar observations are used to compute time correction factors for the 
SeaWiFS bands, which are incorporated into the calibration table. The response of band 7 has decreased by 
3.4%, the response of band 8 has decreased by 9.9%, and the response of bands 1-6 have decreased by less 
than 1%, over the course of the mission. The calibration table will be updated as required by the on-orbit 
performance of the instrument. 


3.1 INTRODUCTION 

The SeaWiFS CVT uses on-orbit calibration data to 
monitor the radiometric stability of the SeaWiFS bands 
over the course of the mission. Monthly lunar are used 
calibrations to track the long-term stability of the radio- 
metric calibration of SeaWiFS. The time series of lunar 
calibrations are used to compute any time correction fac- 
tors required to maintain a stable radiometric sensitivity 
for each of the SeaWiFS bands over the course of the mis- 
sion. The lunar data analysis techniques and preliminary 
results were previously discussed by Barnes et al. (1998 and 
1999a) and by Barnes and McClain (1999). This chapter 
describes how these techniques are used to generate the 
time corrections for the SeaWiFS calibration table. 

The time correction factors are defined in the SeaWiFS 
level -lb calibration equation (discussed in detail in John- 
son et ai. 1999): 

L s ( A) — (C out (A) - C dar k(A)) Ki(g t d, A) 

x (1 + K 2 ( A) (T - T ref )) K 3 ( pxl, A) 
x M(ms, A) q(A) (/3(A) + 7(A) (t — to) 

+ 6(\)(t-t 0 ) 2 ), 

where: 

A is the wavelength; 

Ls is the calibrated at-sensor radiance; 

C ou t Is the counts from sensor output data; 

Cdark Is the dark count from sensor output data; 

K\ is the counts-to-radiance conversion factor (calibra- 
tion coefficient); 


g is the gain; 
d is the detector; 

K 2 is the detector temperature-dependent correction 
factor; 

T is the detector temperature from sensor output data; 
T re f is the reference temperature (20° C); 

Af 3 is the scan modulation correction; 
pxl is the pixel number along scan line; 

M is the half- angle mirror side correction factor; 
ms is the mirror side (a or b); 
a is the vicarious gain; 

P is the constant term in temporal correction; 

7 is the linear (in time) term in temporal correction; 
6 is the quadratic (in time) term in temporal correc- 
tion; 

t is the time tag of sensor output data; and 
to is the reference time for temporal correction. 

The reference time for the temporal correction is the 
time tag of the first SeaWiFS on-orbit image, which was 
obtained on 4 September 1997 at 162630 UTC. The cur- 
rent design of the SeaWiFS level -lb algorithm allows for 
a multisegment quadratic time correction. 

3.2 LUNAR CALIBRATIONS 

Lunar calibrations are performed once per month when 
the moon is at a phase angle (tf) of approximately 7°. This 
phase angle is chosen to maximize the illuminated surface 
of the moon while minimizing the opposition effect, i.e., the 
increase in brightness of sunlight diffusely reflected from a 
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particulate surface near zero phase. Operational consider- 
ations, such as a conflict of the lunar measurement with 
a midnight data downlink, will require that the measure- 
ments be moved on occasion to different phase angles or 
from before full phase to after full phase. 

SeaWiFS operates in a sun-synchronous orbit, cross- 
ing the equator from north to south at local noon. In 
normal operation, the spacecraft is maintained in a nadir 
orientation, using pitch- axis momentum wheels for atti- 
tude control, with a spacecraft pitch rate of 360° per orbit 
(about 0.06° per second). For lunar measurements, the ro- 
tation rate of the momentum wheels is increased and the 
spacecraft is pitched in the opposite direction at a rate of 
approximately 0.15° per second. The maneuver is started 
after the spacecraft crosses the South Pole and is timed so 
that SeaWiFS will view the moon as the spacecraft ground 
track crosses the sublunar point. At the end of the maneu- 
ver, when the spacecraft again points toward the Earth, 
the pitch rate is returned to normal. During the maneu- 
ver, SeaWiFS is oriented so that it scans across the lunar 
surface from west to east in celestial coordinates. 

Because the moon appears to be a stationary object 
during the SeaWiFS measurements, the number of scan 
lines in lunar calibration depends on the pitch rate of the 
instrument and the apparent size of the moon. The pitch 
rate causes SeaWiFS to oversample the moon, resulting in 
approximately 26 scan lines of the moon in a typical lunar 
image that has a diameter of approximately 7 pixels. 


For the calibration data analysis, disk-integrated spec- 
tral radiances are computed for each band from the lu- 
nar images. Prelaunch modeling of simulated lunar im- 
ages showed that disk-integrated spectral radiances pro- 
duce more consistent results than those computed for one 
or a few pixels from the center of the lunar image (Wood- 
ward et al. 1993). In this analysis, the disk-integrated 
radiances are summed over all pixels in the lunar image 
whose brightness is greater than 1% of the peak brightness 
in the image. 

The time series of disk- integrated radiances for the first 
27 lunar calibrations (spanning 798 days from 14 Novem- 
ber 1997 through 21 January 2000) are plotted in Fig. 5. 
The radiances in each band have been normalized by the 
value at the first calibration to show the relative changes in 
the instrument response as functions of time. Much of the 
variation in these time series is due to the changing observ- 
ing geometry of the individual calibrations. Normalization 
of the calibration data for the observing geometry will be 
discussed in the next section. 

3.3 NORMALIZING FACTORS 

Although the surface of the moon remains unchanged 
over the lifetime of SeaWiFS, the radiance from the moon 
varies with the geometry of the observations. As a re- 
sult, the disk-integrated lunar radiances are normalized to 
a common viewing geometry for incorporation into a lu- 
nar calibration time series. These normalizing factors are 
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Table 5. Lunar calibration observing geometry. The notation E3 indicates the calibration was after zero phase; 
[B] indicates before zero phase. The symbol At 0 is the time (in days) since the first image. 


Calibration 

Year 

Date 

At 0 

■Dsm [AU] 

Dim(^m) 

?? [degrees] 

Nm 

1 

1997 

14 Nov. El 

71.26 

0.991602 

0.939681 

6.75 

25.63 

2 


14 Dec. [A] 

100.83 

0.986812 

0.967318 

7.03 

25.35 

3 

1998 

13 Jan. [Al 

130.39 

0.986119 

0.996159 

5.45 

24.57 

4 


10 Feb. E 

159.19 

0.989545 

1.01474 

6.65 

24.29 

5 


12 Mar. E 

188.89 

0.996419 

1.03299 

6.72 

23.62 

6 


12 Apr. [A] 

219.75 

1.00516 

1.03719 

6.66 

24.14 

7 


11 May E 

249.38 

1.01283 

1.02474 

7.10 

25.75 

8 


10 Jun. [a] 

278.87 

1.01790 

1.00304 

6.43 

26.19 

9 


10 Jul. El 

308.36 

1.01919 

0.975564 

5.70 

26.69 

10 


5 Sep. [S 

366.31 

1.01050 

0.932589 

6.52 

27.81 i 

11 


5 Oct. E 

395.73 

1.00238 

0.915500 

6.69 

28.47 

12 


4 Nov. [A] 

425.84 

0.994060 

0.910064 

6.55 

28.16 

13 


4 Dec. El 

455.33 

0.988037 

0.920501 

7.03 

27.24 

14 

1999 

2 Jan. El 

484.89 

0.985745 

0.942602 

6.73 

28.29 

15 


1 Feb. [A] 

514.38 

0.987825 

0.970323 

4.88 

26.79 

16 


2 Mar. El 

544.20 

0.993771 

1.00182 

7.38 

26.34 

17 


31 Mar. OS 

572.73 

1.00158 

1.01719 

7.01 

25.65 

18 


1 May El 

603.38 

1.01015 

1.03724 

6.92 

25.62 

19 


30 May El 

633.21 

1.01647 

1.03694 

7.95 

25.46 

20 


29 Jun. El 

662.83 

1.01928 

1.02351 

7.25 

25.94 

21 


27 Jul. Q2 

691.21 

1.01814 

1.01103 

6.82 

26.15 

22 


26 Aug. [b] 

720.76 

1.01317 

0.984094 

6.72 

26.81 

23 


25 Sep. [B] 

750.32 

1.00553 

0.954810 

6.83 

26.78 

24 


24 Oct. [E 

776.21 

0.997086 

0.929820 

7.25 

27.78 

25 


22 Nov. [B] 

809.28 

0.989973 

0.913459 

6.66 

27.78 

26 


23 Dec. El 

839.72 

0.985987 

0.912118 

9.83 

28.80 

24 


24 Oct. Q2 

776.21 

0.997086 

0.929820 

7.25 

27.78 

25 


22 Nov. [B] 

809.28 

0.989973 

0.913459 

6.66 

27.78 

26 


23 Dec. El 

839.72 

0.985987 

0.912118 

9.83 

28.80 

27 

2000 

21 Jan. El 

869.13 

0.986489 

0.924960 

8.61 

29.02 


based, in large part, on the positions of the spacecraft, 
Earth, sun, and moon computed by the SeaWiFS naviga- 
tion algorithms. The observing geometries for the lunar 
calibrations considered in this analysis are provided in Ta- 
ble 5. Using these values, five normalizing factors are cal- 
culated. These normalizing factors are discussed in detail 
in Barnes et al. (1998 and 1999a). The implementation of 
these normalizing factors is discussed here. 

The first normalizing factor (N\) corrects to a common 
sun-moon distance and varies as l/R 2 , where R is the dis- 
tance between two bodies. The sun-moon distance ( Dsm ) 
is normalized to 1 Astronomical Unit (AUf): 


N i 


Dsm 2 
AU 


t 1 AU = 149.59787066 x 10 6 km. 


(2) 


The second normalizing factor ( N 2 ) corrects to a com- 
mon SeaWiFS-moon distance and varies as l/R 2 . The 
SeaWiFS-moon distance (Dim) is normalized to the mean 
radius of the lunar orbit (Rm = 384401 km): 


N 2 = 


Dim 2 
Rm . 


(3) 


The third normalizing factor (N 3 ) is the illuminated 
fraction of the lunar surface (f\) as a function of phase 
angle. This factor is a linear function of the phase angle, 
with the lunar surface fully illuminated at 0° phase, half 
illuminated at 90° phase, and dark at 180° phase. This 
function is given as: 


/i(tf) = a o + (4) 


where a 0 = 1.0 and ai = -1/180. Because the nominal 
phase angle of the lunar calibrations is 7°, the illuminated 
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area of the lunar surface 
area at 7° : 

N 3 = 


is normalized to the illuminated 


/i(7°) 

/iW 

173.0° 

180.0° - i?' 


(5) 


The fourth normalizing factor (N 4 ) compensates for the 
oversampling of the lunar image during the calibration. As 
discussed previously, the oversampling is a function of the 
pitch rate and the apparent size of the moon. Because the 
spacecraft does not have the use of its horizon sensors dur- 
ing the lunar pitch maneuver, there is increased noise in 
the pitch rate calculated from the attitude control system 
during the maneuver. As a result, the pitch rate is deter- 
mined from the number of scan lines in the lunar image. 
The number of scan lines is computed from the maximum 
extent between the points whose brightness is 1% of the 
peak brightness in the image. For a given calibration, these 
distances are averaged over the eight bands. Because the 
number of scan lines in the lunar calibrations range from 
23-29, these numbers are normalized to a common value 
of 25 scan lines. To compensate for the variation in the 
apparent size of the moon in the images, the number of 
scan lines is normalized to a common SeaWiFS-moon dis- 
tance of the mean radius of the lunar orbit. The resulting 
normalizing factor is: 


N = 25.0/TVm 


( 6 ) 


where N u is the mean number of scan lines in the lunar 
image. 

The fifth normalizing factor (N$) corrects for changes 
in the brightness or reflectance of the moon with phase an- 
gle. The moon has a non-uniform particulate surface with 
large-scale regional variations in reflectance. The nonlam- 
bertian change in the overall reflectance of the lunar sur- 
face with phase angle can be approximated by Hapke’s 
bidirectional reflectance equation (Hapke 1986). Helfen- 
stein and Veverka (1987) used Hapke 's equation and a set 
of six empirically-derived constants, to provide a curve 
of disk-integrated reflectance versus phase angle for the 
moon. This curve is plotted in Fig. 6. A quadratic func- 
tion has been fit to this curve to provide an interpolation 
between the data points (/ 2 ). This interpolation scheme 
is limited to phase angles between 4° and 10°, using the 
function: 

AW = bo + M + M 2 , (7) 

where b 0 = 1.2872531 x 1CT 1 , b x = -6.7007694 x 10" 3 , 
and 62 = 2.1625472 x 10" 4 . The normalizing factor is 
computed relative to the value at 7°, the nominal phase 
angle of the lunar calibrations: 


^5 


/a(7°) 

/aW 

0.092414408 
6q 4- b\ t? + 6 2 1? 2 1 


(8) 


The SeaWiFS lunar calibrations to date have occurred at 
phase angles of 4. 8-9.8°, so the normalizing factor has been 
applied over a narrow range of phase. There are indica- 
tions that the variation in lunar reflectance with phase an- 
gle has a wavelength dependence, which will be discussed 
in subsequent sections. There is also evidence that the 
moon is brighter before full phase than after (Kieffer and 
Anderson 1998). If this evidence is verified, the brightness 
asymmetry would have an effect of 0. 5-1.0% on the value 
of Ns- Approximately 40% of the lunar measurements to- 
date were made before full phase. 

The overall normalizing factor for each lunar measure- 
ment is the product of the five individual factors. This 
multiplicative factor is applied to the disk-integrated lu- 
nar radiances for each of the eight SeaWiFS bands. For 
the lunar calibrations to-date, the value of the overall nor- 
malization factor has ranged from 0.783-1.10, with a mean 
value of 0.924. The time series of radiances normalized to 
the common viewing geometry are plotted in Fig. 7. 

The phase angle is the most important of the geomet- 
ric angular parameters for SeaWiFS lunar measurements. 
The variation of the integrated lunar radiance with phase 
angle is much stronger than any variation with li brat ion 
angle. For libration changes, the loss of visible lunar sur- 
face from one side of the moon is balanced by the gain of 
visible surface on the other side. For disk-integrated spec- 
tral radiances over the course of several months to several 
years, libration is not expected to have a major effect on 
the slope of the time series, but is expected to increase the 
scatter in the data. The overall contribution of libration 
to the SeaWiFS lunar time series remains unknown and is 
not corrected in this analysis. 

3.4 PHASE ANGLE CORRECTIONS 

Two trends are apparent in the plots shown in Fig. 7. 
The first is a systematic variation in the data from one cal- 
ibration to the next that probably arises from an incom- 
plete normalization to a common viewing geometry rather 
than from instrumental effects. The fact that the data for 
bands 1-6 track each other supports this reasoning. The 
second trend is a decrease in the radiometric response of 
bands 7 and 8 with time. 

The time series can be corrected for the incomplete 
normalization over viewing geometry by a second normal- 
ization based on a subset of the SeaWiFS bands. Barnes 
et al. (1998) based this second normalization on the values 
for band 5, while Barnes et al. (1999a) based the normal- 
ization on the average of the values for bands 1-6. After 23 
lunar calibrations, however, the consistency among these 
bands is less than it was for the data sets in Barnes et al. 
(1999a). Barnes and McClain (1999) show that bands 3 
and 4 appear to be changing the least among these bands, 
so the mean value of these bands at each calibration is 
used as the second normalization of the calibration time 
series. The radiances with the second normalization ap- 
plied are plotted in Fig. 8. Identical results are obtained 
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Fig. 8. Disk-integrated lunar radiances, normalized to the first calibration and the mean of bands 3 and 4. 


if the disk-integrated radiances are normalized to the first 
observation in each band and to the mean for bands 3 and 
4 for each calibration. Examining Fig. 8 shows that the 
radiometric responses for bands 1-6 have changed slightly 
over the course of the mission, while the response for band 

7 is down approximately 2.5% and the response for band 

8 is down approximately 8%. It should be noted that cali- 
brations 3, 9, and 15 occurred at phase angles of less than 
6° and that calibrations 19, 26, and 27 occurred at phase 
angles of 8° or more. 

The change in lunar reflectance with phase angle from 
Helfenstein and Veverka (1987) is monochromatic. The 
measurements used as a basis for their lunar reflectance 
model were made at wavelengths from 360-1 ,060 nm (Lane 
and Irvine 1973). Helfenstein and Veverka (1987) used 
the average of those measurements (over wavelength) to 
create a single, best-fit lunar reflectance curve at an un- 
defined wavelength, presumably near 500 nm. The devia- 
tions of the low and high phase angle calibrations shown in 
Fig. 8 appear to vary in magnitude as a function of wave- 
length, indicating that wavelength-dependent phase angle 
effects are still present in the data. Barnes and McClain 
(1999) derived an empirical wavelength-dependent phase 
angle correction from the lunar calibration data that min- 
imizes these deviations. Barnes and McClain (1999) only 
used the data points for calibrations 3, 9, and 15 to com- 
pute the corrections, while the implementation described 
here uses all 27 calibrations to compute the corrections. 

The method for calculating the correction is as follows. 
For each band, a regression line is calculated from the nor- 


malized calibration data to determine the change in the 
output of the band with time. The low and high phase an- 
gle points are excluded from the calculations. Two piece- 
wise fits are computed for each of bands 1-6: one quadratic 
and one linear. Three piecewise fits are computed for 
bands 7 and 8: two quadratic and one linear. For each 
band, the fractional differences between the measured cal- 
ibration data points and the computed trend lines (fa) are 
computed as functions of the phase angles of the measure- 
ments: 

/jM _ £-<*■«)( (9 ) 

where L m is the measured radiance and L c is the com- 
puted radiance. These differences are plotted in Fig. 9, 
along with linear fits to the differences. The slopes of 
these fits, Ci(A), are the wavelength-dependent phase an- 
gle correction factors. Examination of the plots in Fig. 9 
shows there is a correlation between the differences and the 
phase angles over the entire range of phase angles. This 
is the reason that the implementation of the wavelength- 
dependent phase angle correction uses all 27 calibrations 
in computing the correction factors. 

The sixth normalizing factor (7V 6 ) for the lunar data is 
computed relative to the nominal phase angle of 7° : 

Ne = 1.0 - ci(A)(tf - 7.0°). (10) 

The correction factors for each band are given in Table 6. 
Values of N& are typically in the range of 0.987-1.010 for 
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the observations. Because these factors are wavelength- 
dependent, these effects cannot be removed from the cal- 
ibration data by normalizing to the mean of bands 3 and 
4. 


Table 6. Wavelength-dependent phase angle cor- 
rection factors. 


Band 

No. 

A 

[nm] 

Correction Factor c\ 

dm 

1 

412 

-0.0015091569 

2 

443 

-0.0011531493 

3 

490 

-0.00011397443 

4 

510 

0.00011441961 

5 

555 

0.0016632741 

6 

670 

0.0033899319 

7 

765 

0.0041000855 

8 

865 

0.0044748836 


Nq is the final normalization applied to the lunar cal- 
ibration data as part of the series of normalizations dis- 
cussed in Sect. 3.3. The lunar data, with this correction 
applied, is then normalized by the mean values for bands 
3 and 4. The resulting time series are plotted in Fig. 10. 
Comparing Fig. 10 with Fig. 8 shows that the wavelength- 
dependent phase angle correction considerably reduces the 
scatter in the time series. The time series for bands 1-6, 
with an expanded vertical scale, are plotted in Fig. 11. 

3.5 LUNAR TIME CORRECTIONS 

The CVT used the lunar calibration time series shown 
in Figs. 10 and 11 to compute time correction factors for 
bands 1, 2, 5, 6, 7, and 8. The group used linear and 
quadratic fits to the lunar measurements to compute the 
changes in response of those bands with time. There is 
no lunar calibration information prior to the first lunar 
calibration, which was obtained on day 71 after the first 
on-orbit SeaWiFS image was obtained. Consequently, the 
data for each band were renormalized from the plots in 
Figs. 10 and 11 so that the intercepts of the fits have val- 
ues of unity at the time of the first image. The fits were 
sampled at the time of the solar calibrations to facilitate 
comparison of the lunar and solar calibration data. 

The renormalized calibration time series for bands 1-6 
are plotted in Fig. 12. For each of these bands, a single 
linear fit yields the best estimate of the change in response 
of the band with time. These fits are also plotted in Fig. 12. 
Bands 3 and 4 do not show any change in response over the 
course of the mission. The responses of the other bands 
are down: band 1 is down by 0.9%, band 2 is down by 
0.5%, band 5 is down by 0.3%, and band 6 is down by 
0.8%. Even though these changes are small, they could 
produce noticeable effects in the water-leaving radiances 
retrieved from the SeaWiFS ocean data, particularly for 
bands 1 and 6. 


The renormalized calibration time series for all eight 
bands are plotted in Fig. 13. The changes in response 
of bands 7 and 8 can be estimated by fits of two piece- 
wise quadratic functions and one piecewise linear function. 
These functions are also plotted in Fig. 13. This figure 
also shows the linear fits for bands 1-6. Examining Fig. 13 
shows that the change in response for bands 7 and 8 is 3.4% 
and 9.9%, respectively, over the course of the mission. 

The time correction factors for bands 1, 2, 5, 6, 7, and 8 
are computed from the inverses of the fits for these bands. 
There is a single correction for bands 1, 2, 5, and 6 and 
three piecewise corrections for bands 7 and 8. The single 
linear corrections for bands 1, 2, 5, and 6 and the piecewise 
linear correction for bands 7 and 8 are used to extrapolate 
the time corrections for these bands in time beyond the 
date of the last lunar calibration. These time correction 
factors provide a stable radiometric response for each of the 
SeaWiFS bands over the course of the mission. These fac- 
tors are incorporated into the SeaWiFS calibration table, 
which is updated as required by the on-orbit performance 
of the instrument. 

3.6 DISCUSSION 

The CVT is continuing to investigate possible causes of 
the change in response of the SeaWiFS bands with time. 
A comparison of the lunar calibration data with the so- 
lar calibration data shows that the changes in response of 
bands 7 and 8 appears in both the lunar and solar data 
sets. Figure 14 shows the ratio of band 7 to band 8 com- 
puted from the solar measurements, as derived by Eplee 
et al. (2000), superimposed on the ratio computed from 
the lunar measurements plotted in Fig. 13. The solar ratio 
was normalized to have the same value as the lunar ratio 
on day 71 — the day of the first lunar calibration. 

The ratio of the functions fitted to the bands 7 and 
8 lunar data is also plotted in Fig. 14. The plots show 
that the lunar and solar data are consistent and that the 
piecewise functions fitted to the lunar data are consistent 
with the observations. The deviation of the solar ratio 
from the lunar ratio prior to the first lunar calibration is 
caused by changes that occurred in the reflectance of the 
solar diffuser early in the mission (Eplee et al. 2000). The 
agreement between the lunar band 7:8 ratio and the solar 
band 7:8 ratio shows that the changes in response for bands 

7 and 8 arise from changes in the instrument itself, and not 
from artifacts in either the lunar or solar data analyses. 

One scenario for the changes in response of bands 7 and 

8 is as follows. The long wavelength limit for the quantum 
efficiency of silicon photodiodes is about 1 m. At these 
wavelengths, near infrared photons penetrate the surface 
of the material to the point that the response of the diode 
to short-term changes in radiant flux is compromised. It is 
generally believed that, over extended periods of exposure, 
these long wavelength photons can also cause a degrada- 
tion of the diode material. Engineers at the Santa Barbara 
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Research Center, the manufacturer of SeaWiFS, call this 
process “annealing.” The smaller change in response for 
bands 5 and 6 may be due to the same process, just at 
shorter wavelengths. 

The changes in the response of bands 1 and 2 require 
another scenario. Eplee et al. (2000) showed that the re- 
flectance of the SeaWiFS solar diffuser is decreasing with 
time, presumably due to photolyzed organic materials (out- 
gassed from the spacecraft) condensing onto the surface of 


the diffuser. These materials preferentially absorb sunlight 
at shorter wavelengths, so the degradation of the diffuser 
reflectance decreases with increasing wavelength. Such yel- 
lowing of the diffuser was expected, based on experience 
with CZCS. The SeaWiFS primary mirror is located be- 
hind an aperture in its rotating telescope housing, so it is 
not exposed to space in the manner of the solar diffuser. It 
is possible that the effects of the yellowing of the primary 
mirror are just now becoming observable. 
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Chapter 4 


Solar Data Analysis for SeaWiFS Calibration 

Robert E. Eplee, Jr., Robert A. Barnes, 
and Frederick S. Patt 
SAIC General Sciences Corporation 
Beltsville } Maryland 

Abstract 


Daily solar calibrations are used to monitor the short-term variability in the radiometric response of SeaWiFS 
in between the monthly lunar calibrations. The time series of solar observations show the effects of both the 
degradation of the reflectance of the solar diffuser and the change in response of the instrument. The ratio of 
the band 7 time series to the band 8 time series, which minimizes the effects of the diffuser degradation, shows 
that the change in response of bands 7 and 8 in the solar data tracks the change observed in the lunar data. 


4.1 INTRODUCTION 

The SeaWiFS CVT uses on-orbit calibration data to 
monitor the radiometric stability of the SeaWiFS bands 
over the course of the mission. Daily solar calibrations are 
used to track the short-term variability in the instrument 
response in between monthly lunar calibrations. Daily 
detector and gain calibrations are also used to monitor 
the stability of the detectors and postdetector electronics 
within SeaWiFS. Preliminary solar data analysis results 
were previously discussed by Barnes et al. (1999a). This 
chapter describes how the solar, detector, and gain calibra- 
tion data plus corresponding analysis techniques are used 
to monitor the performance of SeaWiFS on orbit. 

4.2 SOLAR CALIBRATIONS 

Solar calibrations are performed once per day when the 
spacecraft is over the South Pole. At this time, the rota- 
tion of the spacecraft causes the sun to rise and set over the 
solar diffuser aperture in the direction of the pitch of the 
spacecraft. The observations record sunlight reflected by 
the uniformly illuminated diffuser. For each solar calibra- 
tion, the mean radiance is computed for ail 8 bands over 
3 pixels and 250 scan lines centered on the image of the 
illuminated diffuser. As the sun changes position in the 
sky over the course of a year, the incidence angle of the 
solar image on the diffuser changes in the direction that is 
perpendicular to the direction of the pitch. This is called 
the “azimuth angle” in the nomenclature of the SeaWiFS 
diffuser (Barnes and Eplee 1996). The azimuthal angles 
for the SeaWiFS solar measurements vary over a range of 
-5° to 6° over the course of a year. The mean solar radi- 
ances are normalized to a common Earth-sun distance of 


one Astronomical Unit and are corrected for the seasonal 
variations in the azimuthal angle of the incident sunlight 
on the diffuser. 

4.3 AZIMUTHAL CORRECTIONS 

Ideally, the bidirectional reflectance distribution func- 
tion (BRDF) of the diffuser would be used to correct for 
the seasonal cycle in the azimuthal angle of the sunlight. 
The BRDF for the SeaWiFS diffuser, however, was only 
determined for bands 4 and 8 under a limited set of view- 
ing angles (Barnes and Eplee 1996). Barnes et al. (1999a) 
fitted the band 8 BRDF data to a second order polynomial 
curve with a value of unity at zero azimuth and a value of 
about 0.95 at 6° on each side of zero to compute an az- 
imuthal angle correction for the solar data. The laboratory 
measurements are incomplete, and the analysis presented 
in this chapter underscores the importance of a detailed 
characterization of the diffuser’s BRDF before launch. 

An alternative approach for correcting the variation in 
the solar azimuthal angle on the diffuser has been devel- 
oped and applied to the solar calibration data. This ap- 
proach assumes that the BRDF of the SeaWiFS diffuser is 
10% of the BRDF of a perfect lambertian diffuser (Barnes 
and Eplee 1996) and incorporates the cosine of the az- 
imuthal angle. The wavelength dependence of the solar 
diffuser (/ 4 ) BRDF is approximated by a rotation in az- 
imuth as a function of wavelength. The resulting correction 
factors are: 

lA 1.0 - (1.0 - cos + (U) 

where: <p is the azimuthal angle of incident sunlight on the 
diffuser; xp is the angle of rotation in azimuth, which is 
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a function of wavelength; and Rb is the ratio of SeaWiFS 
BRDF to lambertian BRDF, which equals 0.1 (Barnes and 
Eplee 1996). The rotation angles for the SeaWiFS bands, 
which have been empirically determined to minimize the 
residual seasonal variations in the solar calibration data, 
are listed in Table 7. 


Table 7. Band-specific rotation angles for azimuth- 
al angle correction. 


A 

Rotation Angie [degrees] 

412 

0.8 

443 

0.7 

490 

0.7 

510 

0.7 

555 

0.6 

670 

0.5 

765 

0.3 

865 

0.0 


With the addition of the wavelength-dependent rota- 
tion angles, the azimuthal angle corrections defined by (11) 
are mathematically equivalent to the polynomial fits em- 
ployed by Barnes et al. (1999a). These adjustments are 
smaller than the 2° spacing of the measurements used to 
create the polynomial curve, and they are based on the as- 
sumption of small alignment differences for measurements 
of the diffuser. 

4.4 CALIBRATION TIME SERIES 

The time series for the first 821 solar calibrations (span- 
ning 864 days from 10 September 1997 through 22 January 
2000) are plotted in Fig. 15. These data have been normal- 
ized to a common Earth-sun distance of 1 AU, but have 
not been corrected for the variation in the solar azimuthal 
angle. The radiances in each band have been normalized 
by the value at the first calibration to show the relative 
changes in the diffuser reflectance as functions of time. The 
time series of solar azimuthal angles, due to the seasonal 
variation in the position of the sun in the sky, is plotted 
in Fig. 16. The azimuthal angle corrections for the eight 
bands, computed using (11), are plotted in Fig. 17. These 
corrections have been applied to the solar calibration data 
plotted in Fig. 18. 

Two trends are apparent in the time series shown in 
Fig. 18. The first is a degradation in the reflectance of the 
diffuser with time as photolyzed organic materials (out- 
gassed from the spacecraft) condense onto the surface of 
the diffuser. These materials preferentially absorb sunlight 
at shorter wavelengths, so the degradation of the diffuser 
reflectance decreases with increasing wavelength. Because 
the correction for the variation in solar azimuthal angle on 
the diffuser is approximate, part of the change in diffuser 
reflectance with wavelength could be a residual BRDF ef- 
fect. The second trend apparent in the data is the changes 


in radiometric responses with time for bands 7 and 8 (dis- 
cussed in detail in Eplee and Barnes 2000). The time series 
do not show any short-term variability in the response of 
the SeaWiFS bands. 

Because the SeaWiFS atmospheric correction algorithm 
uses the ratio of radiances in band 7 to the radiances in 
band 8 for computing the aerosol radiance, and because the 
degradation in the diffuser reflectance decreases with wave- 
length, the ratios of the solar data for band pairs have been 
plotted (Fig. 19). These plots have also been normalized to 
the value of unity at the first calibration. Comparing the 
band pair ratios implies that if there were no degradation 
in the response of bands 7 and 8, the band 7 8 ratio should 
approach unity. The trend in the band 7:8 ratio evident in 
the plot arises from the change in the response of the two 
bands. 

The solar and lunar calibration time series for bands 7 
and 8 are plotted simultaneously in Fig. 20. Comparing 
the solar and lunar data for the two bands provides an 
estimate of the degradation of the diffuser reflectance that 
is independent of wavelength. 

Figure 21 shows the ratios of solar data for bands 7 
and 8 superimposed on the ratio from the lunar calibra- 
tion data (Eplee and Barnes 2000). The solar ratio has 
been normalized to have the same value as the lunar ratio 
on day 71, the day of the first lunar calibration. The plots 
show that the solar and lunar data are consistent and that 
the degradation in the responses of bands 7 and 8 is not 
an artifact of the solar diffuser or of the degradation in the 
diffuser reflectance. The plots also show that the correc- 
tions for the variation in the azimuthal angle (/ 4 ) for bands 
7 and 8 yields data that are consistent with observations 
which are independent of the solar diffuser. 

Time correction factors have been computed for bands 
1, 2, 5, 6, 7, and 8 from the lunar calibration data (Eplee 
and Barnes 2000). These correction factors are applied to 
the solar calibration data plotted in Fig. 22 and to the 
solar calibration band ratios plotted in Fig. 23. The time 
series in Fig. 23 shows that the time corrections for bands 
7 and 8 are stable to within 0.2% over the course of the 
mission. 

4.5 DETECTOR CALIBRATIONS 

Detector calibrations are performed once per day, typi- 
cally two orbits after the solar calibration, using a variation 
of the solar calibration procedure. During the detector cal- 
ibrations, observations of the illuminated solar diffuser are 
collected individually for the four detectors in each band. 
Figure 24 shows the time series of the detector calibrations 
for band 7, and Fig. 25 shows the time series of detector 
calibrations for band 8. These plots show that the changes 
in response of the individual detectors in bands 7 and 8 
are consistent and track the drifts of the bands as a whole. 
The changes in response of bands 7 and 8 are not due to 
changes in individual detectors within either band. 
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Fig. 15. Solar calibration diffuser reflectances. 



Fig. 16. Azimuthal angles of the image of the sun on the solar diffuser. 
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Table 8. SeaWiFS prelaunch gains relative to Gain 1. 


Band No. 

Detector No. 

Gain 1 

Gain 2 

Gain 3 

Gain 4 

1 

1 

1.00000 

1.01562 

1.00506 

1.01046 


2 

1.00000 

1.98647 

1.32275 

1.68180 


3 

1.00000 

1.98621 

1.32010 

1.68059 


4 

1.00000 

1.99245 

1.31944 

1.67965 

2 

i 

1.00000 

1.99157 

1.31919 

1.68193 


2 

1.00000 

1.98408 

1.31888 

1.68100 


3 

1.00000 

1.98892 

1.31722 

1.68275 


4 

1.00000 

1.02419 

1.00806 

1.01613 

3 

1 

1.00000 

1.05087 

0.999861 

1.03091 


2 

1.00000 

1.98867 

0.894868 

1.68113 


3 

1.00000 

1.98857 

0.896848 

1.68113 


4 

1.00000 

1.99015 

0.895656 

1.67992 

4 

1 

1.00000 

1.99007 

0.786856 

1.68143 


2 

1.00000 

1.98781 

0.788264 

1.68176 


3 

1.00000 

1.98821 

0.788229 

1.68360 


4 

1.00000 

1.04348 

0.985507 

1.02899 

5 

1 

1.00000 

1.04412 

0.985294 

1.02941 


2 

1.00000 

1.99080 

0.641401 

1.59637 


3 

1.00000 

1.98940 

0.643069 

1.59636 


4 

1.00000 

1.98948 

0.641937 

1.59425 

6 

1 

1.00000 

1.99099 

0.362975 

0.662879 


2 

1.00000 

1.98762 

0.364917 

0.665017 


3 

1.00000 

1.99050 

0.364214 

0.665150 


4 

1.00000 

1.03333 

0.966673 

0.979392 

7 

1 

1.00000 

1.04918 

0.966610 

0.968238 


2 

1.00000 

1.98794 

0.310052 

0.575541 


3 

1.00000 

1.98864 

0.312990 

0.575752 


4 

1.00000 

1.99055 

0.310930 

0.576902 

8 

1 

1.00000 

1.99170 

0.260256 

0.498881 


2 

1,00000 

1.98912 

0.261251 

0.499517 


3 

1.00000 

1.99130 

0.259961 

0.498223 


4 

1.00000 

1.04110 

0.958904 

0.972603 


Table 9. SeaWiFS on-orbit gains relative to Gain 1. 


Band No. 

Detector No. 

Gain 1 

Gain 2 

Gain 3 

Gain 4 

1 

1 

1.00000 

1.01384 

1.00474 

1.01070 


2 

1.00000 

1.98483 

1.32261 

1.68071 


3 

1.00000 

1.98667 

1.31367 

1.67983 


4 

1.00000 

1.99250 

1.31748 

1.67960 

2 

1 

1.00000 

1.99404 

1.31923 

1.67573 


2 

1.00000 

1.98330 

1.31885 

1.68066 


3 

1.00000 

1.99005 

1.31212 

1.68243 


4 

1.00000 

1.02201 

1.00803 

1.01613 
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Table 9. (cont.) SeaWiFS on-orbit gains relative to Gain 1. 


Band No. 

Detector No. 

Gain I 

Gain 2 

Gain 3 

Gain 4 

3 

1 

1.00000 

1.04317 

0.999249 

1.02541 


2 

1.00000 

1.96727 

0.895505 

1.68115 


3 

1.00000 

1.99016 

0.893043 

1.68186 


4 

1.00000 

1.99132 

0.896420 

1.68058 

4 

1 

1.00000 

1.99264 

0.787118 

1.67559 


2 

1.00000 

1.98548 

0.787691 

1.68128 


3 

1.00000 

1.98929 

0.785179 

1.68326 


4 

1.00000 

1.04064 

0.985578 

1.02899 

5 

1 

1.00000 

1.03909 

0.984992 

1.02914 


2 

1.00000 

1.98856 

0.642293 

1.59672 


3 

1.00000 

1.99075 

0.640013 

1.59658 


4 

1.00000 

1.99070 

0.643690 

1.59483 

6 

1 

1.00000 

1.99330 

0.363494 

0.666026 


2 

1.00000 

1.98551 

0.365785 

0.668350 


3 

1.00000 

1.99166 

0.362301 

0.665366 


4 

1.00000 

1.03082 

0.966881 

0.980609 

7 

1 

1.00000 

1.04439 

0.955551 

0.967513 


2 

1.00000 

1.98649 

0.311277 

0.579590 


3 

1.00000 

1.98952 

0.311649 

0.575681 


4 

1.00000 

1.99139 

0.314473 

0.576669 

8 

1 

1.00000 

1.99373 

0.260326 

0.503102 


2 

1.00000 

1.98696 

0.262634 

0.503937 


3 

1.00000 

1.99227 

0.258170 

0.498930 


4 

1.00000 

1.03832 

0.959105 

0.972930 



Fig. 25. Detector calibration time series for band 8. 
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4.6 GAIN CALIBRATIONS 

Gain calibrations are performed once per day, typically 
immediately following the solar calibrations, and are used 
to track the stability of the four gains of each detector 
in all eight bands. During the calibration, a square-wave 
electronic pulse is applied to the post detector electronics 
during a portion of the SeaWiFS scan when the detectors 
are not illuminated. The gain of the detectors is varied over 
the duration of the pulse. The ratios of the output from the 
detectors at each gain, relative to the output for the unity 
gain, are the gain ratios. The gain ratios, computed from 
the on-orbit calibration pulse data, have been stable over 
the course of the mission and are consistent with the ratios 
measured during the prelaunch instrument recalibration 
(Johnson et al. 1999). The prelaunch ratios are shown in 
Table 8 and the on-orbit ratios are shown in Table 9. A 
comparison of the two shows the gain ratios are stable to 
within 0.1%, which eliminates the postdetector electronics 
as a cause for changes in the response of bands 7 and 8. 


4.7 DISCUSSION 

Despite the degradation of the solar diffuser reflectance 
with time, the SeaWiFS solar calibration data show the 
same changes in radiometric responses of bands 7 and 8 
that were first observed in the lunar calibration data. This 
change in response is not caused by any artifacts in the so- 
lar data analysis. The solar data show no short-term varia- 
tions in the response of the SeaWiFS bands between lunar 
calibrations. Because of the degradation of the diffuser re- 
flectance, the small drifts in the response of bands 1, 2, 
5, and 6, which were observed in the lunar data, are not 
visible in the solar data. 

The analyses of the on-orbit detector and gain cali- 
bration data imply that the source of the changes in the 
radiometric response of bands 7 and 8 is found either prior 
to the detectors in the optical path of the instrument, or 
it is the result of a change that occurs simultaneously to 
all four detectors in each band. 
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Chapter 5 


Vicarious Calibration of SeaWiFS Band 7 

Wayne D. Robinson 
SAIC General Sciences Corporation 
Beltsmlle, Maryland 

Menghua Wang 

University of Maryland Baltimore County 
Baltimore, Maryland 

Abstract 

A technique for vicariously calibrating the gain in the SeaWiFS band 7 is described in this chapter, as well as 
the details of its implementation. The technique has been found to provide a consistent estimate of the gain 
over the life of the SeaWiFS Project. 


5.1 INTRODUCTION 

This chapter describes the current method used to vi- 
cariously determine the gain factor to apply to the band 
7 SeaWiFS radiances. Because of their negligible water- 
leaving radiance components, band 7, at 765 nm, is used in 
combination with band 8 at 865 nm to determine the aero- 
sol optical properties of the air column (Gordon and Wang 
1994b), from which the aerosol radiances in the remainder 
of the SeaWiFS bands can be deduced. Determining the 
band 7 gain properly is important so that all the other 
bands can be calibrated and used properly. 

The problem of vicarious calibration of ocean color sen- 
sors was discussed by Gordon (1998), who proposed the 
following requirements for achieving a good calibration: 

1. The calibration should be made in a cloud-free air 
mass with a maritime aerosol having an optical 
thickness of less than 0.1, and 

2. The water-leaving radiances must be uniform over 
the area in question. 

In the case of band 7 (and 8), these requirements should be 
easily met, as long as highly turbid or high chlorophyll con- 
centration waters are avoided and clear areas are selected. 
Many of the quantities that need to be measured to prop- 
erly perform the calibration are available in the form of 
SeaWiFS radiances and ancillary data. Specific knowledge 
of the aerosol type is not available, but this uncertainty is 
reduced by using a region where the aerosol properties are 
well known. 

Section 5.2 investigates the theory behind the calibra- 
tion of band 7, Sect. 5.3 examines the validity of assuming 
band 8 has the correct calibration, and Sect. 5.4 discusses 


the technique for evaluating the gain in more detail, using 
the measurements available to the SeaWiFS Project. The 
results and conclusions of these analyses are presented in 
Sects. 5.5 and 5.6. 


5.2 THEORY 

Gordon and Wang (1994a) derived a value, €, based on 
the single scattering aerosol reflectance, p as : 


e(A, 865) 


Pas(A) 

* 7 ,( 865 )’ 


( 12 ) 


which is a constant for any one type of aerosol atmos- 
phere. Computations of e have been made for a variety 
of aerosol types including oceanic, maritime, coastal, and 
tropospheric, and at a variety of relative humidities: 50, 
70, 90, and 99%. SeaWiFS bands 7 and 8 were chosen 
so that the e(765,865) could be computed. The equation 
for upwelling radiance through a clear atmosphere over the 
ocean at the top of the atmosphere (TOA) is: 


Lt{X) — L r { A) -f L a (A) T L ra ( A) 

+ T(X)L 9 (X) T t(X) Lf(X) (13) 
+ t(X)L w (X), 

where L t ( A) is the total upwelling radiance, L r ( A) is the 
Rayleigh radiance, L a (A) is the radiance arising from aero- 
sol scattering, L ra (A) is the radiance arising from the in- 
teraction of molecular and aerosol scattering, T(X)L g (X) is 
the glint radiance arising from the specular reflection of the 
sun off the water surface, *(A)L/(A) is the radiance from 
foam (or whitecaps), and t(X)L w (X) is the water-leaving 
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radiance. Note that T( A) is the direct transmittance and 
t(X) is the diffuse transmittance of the atmosphere. L r (A) 
can be determined accurately with a knowledge of the sur- 
face pressure (Gordon et al. 1988). Areas where the sun 
glint is significant can be predicted and avoided by viewing 
away from the point of specular reflection. The whitecap 
radiance can be well estimated at low wind speeds (Gor- 
don and Wang 1994b) and can be avoided at higher wind 
speeds as the wind field is an available product. Lw(X) at 
the 765 and 865 nm bands of SeaWiFS can be considered 
to be zero in nonturbid, low chlorophyll waters; thus, (13) 
can be simplified for bands 7 and 8 to: 

L t (X) - L r (X) - t(X)Lf(X) = L a (X) + L ra (X). (14) 

Gordon and Wang (1994a) determined a relationship be- 
tween the L 0 (A) 4- L ra { X) terms and the single scattering 
radiance, L oa (A), for the aforementioned aerosol models; 
thus, because L r (A) and t(X)Lj(X) are known, the total 
radiance for bands 7 and 8 can be used in (12) to derive e. 

Over open ocean areas, it can be assumed that a mar- 
itime aerosol is usually present and, thus, such a site has a 
known e value. For sites like this, it is then possible to de- 
termine the gain in band 7 through a process of adjusting 
the gain until (12) produces a maritime e value. 

5.3 BAND 8 ACCURACY 

The basic method described in Sect. 5.2 can be applied 
to SeaWiFS observations to determine the gain for band 
7. Before this method is described, a tacit assumption is 
made that the gain in band 8 is already correct, but if it 
is not, the band 7 gain can absorb any small offsets in the 
band 8 gain. In fact, band 8 does have a laboratory cal- 
ibration (Johnson et al. 1999) and should be reasonably 
close to the correct value (Barnes et al. 1999b); however, 
the calibration can change as a result of the stresses asso- 
ciated with launch. Also, the calibration being performed 
is not strictly an absolute calibration, but a vicarious cali- 
bration which includes a calibration of the sensor and the 
processing algorithms as a whole. Studies of the errors in- 
troduced by gain errors (Gordon 1998, Wang and Franz 
2000) indicate that this assumption introduces only minor 
errors. An analysis of the error introduced in the band 7 
calibration due to an incorrect band 8 calibration is pre- 
sented here. 

If the single-scattering approximation is used, the e 
value used in the band 7 gain calculation can be expressed 
as: 


[L t (765)G 7 -L r (765)]F 0 (865) 
e(765, 865) « j Lt ( 865 ) G8 _ L r (865)j F 0 (765)’ 


(15) 


where G 7 and Gg are the gain factors for bands 7 and 8, 
and F 0 (X) is the solar irradiance. Assuming that the band 


7 gain can compensate for any errors in the band 8 gain, 
then (15) can be rewritten as 


e(765, 865) 


[L t (765)G 7:8 -L r (765)]F 0 (865) 

[L f (865) - L r (865)] F 0 (765) ’ J 


where G 7 :s is the combined gain used as the band 7 gain 
in this analysis. If L r (A) is small, then 



(17) 


and the assumption would work perfectly; however, L r (A) 
is significant relative to Lt( A). 

By converting (15) into the form of (16), the actual 
value of Gr.s is seen to be: 


G 7 : 8 = G7h+!j4^(l-h), (18) 

where h is a shorthand notation for 

= L t ( 865) - L r (865) 

Lt(865) Gg - L r (865) V } 

The actual value of Gr.g is not of interest, but the 
amount of variation, specifically, the standard deviation in 
Grs, i.e., cr(Gr.g), for the normal range of SeaWiFS data is 
important. This would indicate the error to be expected in 
the gain applied to band 7 for an average SeaWiFS pixel. 

Gr.$ was evaluated for a typical G AC pass of data using 
a nominal band 7 gain (G 7 ) of 0.95, and values of Gg rang- 
ing from 0.95-1.05 , i.e., band 8 gain changes of from -5% 
to -1-5%. For the 11,778 pixels that could be processed to 
get L W n data, the cr(G 7 :e) is shown in Table 10. 

Table 10. Error in the band 7 gain. Column 1 
is the actual gain in band 8 for which a gain of 1 
is assumed. Column 2 is the error in the band 7 
gain expected (i.e., an extra source of noise) and is 
a result of assuming the band 8 gain is 1. 


^8 

^(G 7 : s) 

0.95 

0.00205 

0.98 

0.00179 

0.99 

0.00091 

1.01 

0.00094 

1.02 

0.00191 

1.05 

0.00489 


The error in the present estimate of band 7 gain is 0.008 
(Sect. 5.5), so unless the error in the band 8 gain is large, 
5% or more, this error is probably acceptable at this time. 
If the gain in band 8 needs to be changed, tests indicate 
that a 3% increase (or decrease) in the band 8 gain requires 
a corresponding band 7 increase (or decrease) of 2.2% to 
get the same e value. 
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Fig. 26 . An image of the Hawaiian Islands showing the location of the MOBY site used to derive the 
band 7 gain. The islands of Hawaii, Maui, Molokai, Lanai, and Oahu are visible in the image with a 
black outline of their coasts. The white box shows the 3x3 pixel area centered on the MOBY location 
at 20.828°N and 157.19°W. 


5.4 CALIBRATION METHOD 

This study was conducted with SeaWiFS data taken 
at the MOBY site located west of the Hawaiian island of 
Lanai (Fig. 26). The site is in the open ocean and has 
a consistent maritime aerosol type. This region was also 
chosen because it has had good coverage by the full, 1 km 
resolution LAC data throughout the mission and it coin- 
cides with the area where the vicarious calibration of the 
other bands is done (Eplee and McClain 2000). 

The choice of the MOBY site allowed the use of 724 
LAC data sets for the study. Such a large number makes 
it possible to be more selective of the atmospheric condi- 
tions in each observation and still retain a large sample of 
observations. The site was chosen to have a square area of 
3x3 LAC pixels, which are averaged together to determine 
the e value for that observation. The t values found in good 
pixels are averaged and then matched against the expected 
maritime c value. The gain in band 7 that produces this e 
value is the gain determined for that observation. 

Although the region for this study was selected to have 
a constant maritime aerosol value, this still leaves four mar- 
itime models to choose from with varying humidity condi- 
tions. If there was a measure of the aerosol value at the 
site, the exact maritime model could be determined. The 
humidity is available in the ancillary data but there may be 


little correspondence between the actual humidity and the 
aerosol model based on that humidity. So, in this study, 
the e value used for the site is the average of the e values 
from the four maritime models. 

The computation of e is performed by running the Sea- 
WiFS operational level -2 program. Some changes are made 
in the standard flagging controls to ensure that only very 
clear sites are used. First, the cloud albedo threshold 
is changed from 1.1 to 0.9 so that more stringent cloud 
screening is used to ensure the selection of cloud-free pixels. 
The following pixel exclusion conditions are also added: a 
satellite zenith angle limit of 56° , a solar zenith angle limit 
of 70° , and a mask of pixels containing excessive stray light. 
Also note that the calibration table used in the processing 
has had the time dependence in bands 7 and 8 removed 
(Eplee and Barnes 2000). The result of the level -2 pro- 
cessing is a 3x3 field of e values and other parameters. 
The unmasked parameter values are averaged to produce 
an observation for that site and data set. The averaged e 
value is used to match the four- model maritime e value. 

Additional screening tests are applied to the observa- 
tions. Observations are considered only if more than five 
of the points remained unmasked. This screen is designed 
so that a relatively large clear region is used for each ob- 
servation and to reduce the effects of noise by averaging 
a number of LAC pixel values of t. Tests were run to see 
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Fig. 27. The gain in band 7 plotted as a function of the number of days since the start of SeaWiFS 
operations on 4 September 1997. The solid line is at the mean gain 7 value and the dashed lines are 1 a 
away from the mean. 



Fig. 28. A histogram of the distribution of band 7 gain at the MOBY site. 


if a requirement of 100% of the points would reduce the 
error in the band 7 gain; it did not, instead however, it 
significantly reduced the number of usable observations. 

In the processing, the effect of whitecap radiance is 
accounted for using a relationship between the wind speed 
and the whitecap radiance (Gordon and Wang 1994a). The 
data used to derive this relationship shows a wide variance 
as the wind speed increases. In addition, during some tests, 
it was found that there may be a weak correlation between 
wind speed and the band 7 gain. Considering this, another 
screen is imposed on the observations that keeps only the 
observations where the wind speed is less than 8ms” 1 . 

The clear air screening was improved using the criteria 
of Gordon (1998) on the aerosol optical thickness. Only 
observations having an aerosol optical thickness less than 


0.1 are used to determine the band 7 gain. A low limit 
of 0.03 is also imposed on all of the observations. A fi- 
nal screen is imposed on the observations to remove any 
statistically bad observations. The standard deviation of 
the band 7 gain is derived for the remaining observations 
and any observations that are more than 2a away from the 
mean are discarded. 

5,5 RESULTS 

The technique described above is applied to available 
SeaWiFS LAC data using the operational level *2 process- 
ing program. Figure 27 is a plot of the band 7 gain de- 
termined for the MOBY site as a function of the num- 
ber of days since SeaWiFS became operational. Out of 
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the 724 LAC data sets covering the Hawaii area, 89 data 
sets satisfy the screening tests; the mean band 7 gain is 
0.940 ±0.008. The error in the band 7 gain of less than 
1% translates into errors in the normalized water-leaving 
radiance in the other bands of well below 10% (Gordon 
1998). This assumes that there is no vicarious calibration 
of the other bands. The vicarious calibration should com- 
pensate for any possible gain 7 errors. Figure 28 shows 
the histograms of the band 7 gain at the MOBY site. 


5.6 CONCLUSIONS 

A technique for calibrating the gain in SeaWiFS band 
7 was derived in accordance with the requirements set out 
by Gordon (1998) and using the resources available to the 
SeaWiFS Project. A reasonably good estimate of the band 
7 gain was made for the current data. This method will be 
applied as new calibration points and algorithm improve- 
ments become available. 
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Chapter 6 


MOBY Data Analysis for the Vicarious 
Calibration of SeaWiFS Bands 1—6 

Robert E. Eplee, Jr. 

SAIC General Sciences Corporation 
Beltsville, Maryland 

Charles R. McClain 
NASA Goddard Space Flight Center 
Greenbelt y Maryland 

Abstract 

The CVT performs the vicarious calibration of SeaWiFS by comparing normalized water-leaving radiances re- 
trieved from SeaWiFS imagery with contemporaneous measurements of normalized water-leaving radiances from 
MOBY. This is a system-level calibration that incorporates the performance of both the sensor and the atmo- 
spheric correction algorithm. The vicarious gains for bands 1-6 depart from unity over a range of approximately 
1-4%, depending on the band, which is consistent with the uncertainty in the prelaunch calibration of SeaWiFS. 


6.1 INTRODUCTION 

Variations in the radiometric response of the eight Sea- 
WiFS bands with time, and uncertainties in the atmo- 
spheric correction algorithm, require a mission-long vicari- 
ous calibration program to monitor the performance of the 
sensor system, instrument plus atmospheric correction al- 
gorithm (Evans and Gordon 1994), to meet the radiometric 
constraints on the ocean color data set (McClain et al. 1992 
and 1998). The temporal stability of the on-orbit calibra- 
tion of the instrument is monitored with a time series of 
lunar calibration data (Eplee and Barnes 2000). The CVT 
is using data from MOBY (Clark et al. 1997), deployed off 
of Lanai, Hawaii, for the vicarious calibration of SeaWiFS 
bands 1-6. 

The CVT performed the vicarious calibration of Sea- 
WiFS by comparing normalized water-leaving radiances, 
Lwn (Gordon and Clark 1981), measured by MOBY with 
Lwn retrieved by SeaWiFS from contemporaneous over- 
flight images of the buoy site. The vicarious calibration 
process adjusts the prelaunch calibration gains to minimize 
the difference between the SeaWiFS and MOBY Lwn in 
bands 1-6. Because no methodology has been developed 
to vicariously calibrate bands 7 and 8, the gains for these 
bands are adjusted to optimize the atmospheric correction 
in the vicinity of MOBY (Robinson and Wang 2000). 

The vicarious calibration is independent of the time 
correction derived from the lunar calibrations (Eplee and 
Barnes 2000). The CVT has derived a set of system gains 
which, when applied to the SeaWiFS calibration, yields 


values for Lwn measured by SeaWiFS and MOBY that 
agree to better than 1%. These gains are defined in a sim- 
plified version of the SeaWiFS level -lb calibration equa- 
tion: 

L s { A) = (C out (A) - C dark (A)) Ki(g,d,\) 

x K 2 { A) or(A) [/3(A) + 7(A) (20) 

x (t - t 0 ) + *(A)(t-t 0 ) 2 ], 

where: 

A is the wavelength of measurement; 

Ls is the calibrated at-sensor radiance; 

C ou t Is the counts from sensor output data; 

Cdark is the median vdL ^ ue of dark count from sensor out- 
put data; 

K\ is the counts to radiance conversion factor (calibra- 
tion coefficient); 
g is the gain; 
d is the detector; 

K 2 is the additional calibration factors; 
a is the vicarious gain; 

0 is the constant term in temporal correction; 

7 is the linear (in time) term in temporal correction; 
6 is the quadratic (in time) term in temporal correc- 
tion; 
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t is the time tag of sensor output data; and 

to is the reference time for temporal correction. 

The reference time for the temporal correction is the time 
tag of the first SeaWiFS on-orbit image, which was ob- 
tained on 4 September 1997 at 162630 UTC. The full level- 
lb calibration equation is presented in Eplee and Barnes 
( 2000 ). 

In some of the SeaWiFS level-la data sets, the dark 
count for particular scan lines may have spuriously high 
values. In addition, the dark count radiance for some 
bands is midway between two digitization levels. These 
two effects can give rise to a striping effect in the level -2 
products. To avoid this striping, the median value of the 
dark counts for each band over each scene is computed and 
subtracted from each scan line in the scene. 

The determination of the temporal correction factors 
is presented in Eplee and Barnes (2000) and discussed in 
Barnes et al. (1999a). The vicarious calibration strategy 
employed by the CVT assumes that the temporal correc- 
tions to the instrument calibration yield stable top-of-the- 
atmosphere radiances from SeaWiFS. This paper discusses 
the current strategy employed by the CVT to determine 
the vicarious gains, a, for each band. 

6.2 NEAR-INFRARED CALIBRATION 

Because open ocean reflectances are low («2%), ap- 
proximately 90% of the top-of-t he- atmosphere signal ob- 
served by SeaWiFS over the oceans is due to Rayleigh 
scattering of sunlight and to aerosol radiance within the 
atmosphere. The SeaWiFS atmospheric correction algo- 
rithm must remove this atmospheric signal to yield the 
water-leaving radiances. 

The atmospheric correction algorithm estimates the 
aerosol radiance L A (which is equal to L a + L ro ), for bands 
7 and 8, and extrapolates L A in the other SeaWiFS bands 
using the ratio of L A in band 7 to that in band 8 (Gor- 
don and Wang 1994a). This ratio is called e [e(765, 865)]. 
Currently, the vicarious gain for band 8 is defined to be 
unity. 

The calibration of band 7 is accomplished by acjjust- 
ing the gain for band 7 so that the c value has the ex- 
pected value for a set of open-ocean scenes in the vicinity 
of MOBY. This procedure is discussed in Robinson and 
Wang (2000). The gain for band 7 derived for the vicari- 
ous calibration discussed here is 0.946. 

6.3 VISIBLE BAND CALIBRATION 

In performing the vicarious calibration, the CVT has 
produced a match-up data set of simultaneous observations 
of Lwn from SeaWiFS and MOBY observations. The cur- 
rent vicarious calibration data set contains 125 matchups 
spanning a 906 day time range from 19 September 1997 
through 13 March 2000. 


6.3.1 SeaWiFS Data Selection 

The SeaWiFS observations are mean water-leaving ra- 
diances ( Lw ) computed for 3x3 pixel regions centered on 
the pixel containing MOBY, where at least five pixels in 
the region pass the exclusion criteria: 

a) Land; 

b) Clouds and ice; 

c) Sun glint; 

d) Stray light; 

e) Total radiance above the knee of the bilinear gain; 

f) Low water-leaving radiance in band 5; 

g) Atmospheric correction algorithm failure; 

h) Scan angle greater than 45°; 

i) Satellite zenith angle greater than 56°; 

j) Solar zenith angle greater than 70°; 

k) Turbid water; 

l) Coccolithophore; and 

m) Aerosol optical depth in band 8 greater than 0.1. 
These criteria are based on standard quality control masks 
and flags, computed on a pixel- by-pixel basis. It should be 
noted that some of these criteria are not directly applicable 
to observations obtained around MOBY, such as turbid 
water and coccolithophores, but are included to maintain 
consistency with other SeaWiFS match-up analyses such as 
Bailey et al. (2000). Sun glint in the SeaWiFS scenes can 
be interpreted by the atmospheric correction algorithm as 
aerosol radiance. To avoid sun glint contamination of the 
match-up data, an upper limit of 0.1 is set on the aerosol 
optical depth in band 8 for valid SeaWiFS retrievals. The 
imposition of this limit results in the loss of several match- 
up scenes during the summer. 

Alternative statistical measures to the mean value 
(which have been considered for determining the optimum 
value of Lw in a scene), include the median value of the 
pixels, the value of the central pixel, and the mode of the 
pixels. The median value can be affected by outliers. Be- 
cause the standard deviation of the mean is typically small, 
the central pixel does not provide a better value than the 
mean. The mode may provide the best estimate of L w , 
but it is difficult to compute for nine pixels. Consequently, 
the mean value is used in the vicarious calibration. 

The mean Lw are converted to Lwn for the match-up 
comparison as discussed in Sect. 6.3.3. The time series of 
mean SeaWiFS Lwn values for the match-up scenes are 
plotted in Fig. 29. The primary source of noise in the 
plots is the variation in the atmospheric correction of the 
SeaWiFS data. 

6.3.2 MOBY Data Selection 

MOBY measures downwelling irradiances and upwelling 
radiances over the wavelength range of 340-900 nm at a 
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subnanometer resolution with two spectrometers coupled 
by a dichroic beamsplitter. The beamsplitter gives the 
blue spectrometer a bandpass of 340-600 nm and the red 
spectrometer a bandpass of 630-900 nm. The potential for 
stray light is greatly reduced by splitting the visible spec- 
trum at the beginning of the water absorption region, be- 
cause most of the short wavelength energy is diverted from 
the entrance slit of the long wavelength spectrometer. The 
splitting also allows the spectrometers to be optimized, in 
terms of free spectral range and integration times, for the 
two distinctive spectral domains. 

The MOBY observations are mean water-leaving ra- 
diances measured over 30 min. intervals centered on local 
noon, the satellite overpass time. In-water measurements 
are made at depths of 2, 5, and 9 m. These measurements 
are used to compute diffuse attenuation coefficients at each 
depth and, thus, to derive Lw at the surface. Surface irra- 
diance ( E s ) measurements are also made. For the MOBY 
spectra to be considered valid, the diffuse attenuation co- 
efficients computed for each depth must be consistent with 
each other. The calibrated MOBY spectra are convolved 
with the SeaWiFS relative spectral response functions for 
use in the match-up analysis. The processed MOBY data 
for a given day includes Lw f° r bands 1-6 and E s for 
bands 1-8. Estimates of E 8 from the subsurface data are 
currently not provided. 

High calibration accuracy for the MOBY data requires 
that the buoy in the water is swapped out for refurbish- 
ment and recalibration approximately every three months. 
To ensure continuous data, the MOBY Project maintains 
three buoys, one in the water and two undergoing refur- 
bishment. The current vicarious calibration uses match- 
up data from seven MOBY deployments. The mean Lw 
is converted to Lwn for the match-up comparison as dis- 
cussed in Sect. 6.3.3. The time series of MOBY Lwn are 
plotted in Fig. 30. The plots show that the MOBY data 
are stable from one deployment to the next. 

Shortly after the launch of SeaWiFS, the above- water 
E s detector on MOBY failed. The E a measurements from 
MOBY for October and November of 1997 are invalid. 
This failure is one of the reasons that the CVT uses com- 
puted E s rather than measured E s for estimating MOBY 
Lwn values in the match-up comparisons of the vicarious 
calibration. 

6.3.3 Match-Up Analysis 

In performing the vicarious calibration, the CVT eval- 
uated matchups of Lwn between SeaWiFS and MOBY. 
Initially, SeaWiFS retrievals of Lwn were compared with 
Lwn computed from MOBY measurements of Lw and E 3 : 

LwnW = ^g-^ Fo(A), (21) 

where Fq is the solar constant and 9o is the solar zenith 
angle. Because of problems with the E s measurements, 


such as the detector failure mentioned above, this approach 
was not used in the match-up analysis. 

For the match-up analysis, normalized water-leaving 
radiances were computed from the SeaWiFS retrievals of 
Lw and from the MOBY measurements of Lw by using 
the atmospheric diffuse transmittance, <(A): 


LwnW 


Lw{ A, #o) 
cos(#o) £(A) 


(22) 


The definition for t(A) is: 


t( A) — exp 


— (0.5r r (A) + t oz (A) + Kg) 
cos (00 ) 


(23) 


where: 

r r is the computed Rayleigh optical thickness; 
r oz is the ozone optical thickness from the SeaWiFS 
ancillary data; and 

K a is the aerosol effects (0.0054) estimated for an aero- 
sol optical thickness of 0.1 and the maritime aerosol 
model for 90% relative humidity. 

The approximation of the aerosol effects with a value in- 
dependent of wavelength is made because La has only a 
weak dependence on wavelength for the maritime aerosol 
model. SeaWiFS retrievals of Lwn can be compared di- 
rectly with the Lwn computed from the MOBY data using 
(22). The vicarious calibration, however, is based on the 
computation of Lwn from Lw using (22) for both Sea- 
WiFS and MOBY to minimize the uncertainties due to 
the approximations made in (22). 

As a check on the Lwn match-up analysis, the com- 
parisons were also performed for Lw- Because the MOBY 
observations are centered on the satellite overpass times, 
the Lw match-ups yielded essentially the same results at 
the Lwn matchups, with equivalent levels of noise in the 
data. 

In performing the vicarious calibration, the vicarious 
gains, a(A) in (20), were adjustd to optimize the agree- 
ment between the Lwn retrieved by SeaWiFS with the 
Lwn measured by MOBY in each band. To avoid any 
residual seasonal variations in solar illumination, the ra- 
tios of the SeaWiFS values to the MOBY values (Lyjfi) 
were computed in each band for each scene. Because the 
distribution of the ratios of Lwn I s more log-normal than 
normal, the vicarious gains were adjusted until the geo- 
metric mean of the ratios for each band was essentially 
unity. 

The results of the vicarious calibration, shown in Ta- 
ble 11 and plotted in Fig. 31 as functions of time, were de- 
rived using an atmospheric correction algorithm which in- 
corporates the near-infrared reflectance correction of Siegel 
et al. (2000). Other estimators of the L$$ over the match- 
up scenes were considered (i.e., the arithmetic mean, the 
median, and the mean of the center quartiles), but the use 
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Fig. 29. SeaWiFS normalized water-leaving radiances for match-up scenes 
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Fig. 30. MOBY normalized water-leaving radiances for match-up scenes. 
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Fig. 31. Vicarious calibration matchups versus time. 
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Table 11. Results of the SeaWiFS vicarious calibration. The MOBY and SeaWiFS radiances are in units of 


mW cm 2 sr 1 1 . 


A 

w 

MOBY 

Radiance 

SeaWiFS 

Radiance 

Mean Ratio 
(SeaWiFS/MOBY) 

Vicarious 

Gain 

412 

1.8263525 

1.8263380 

0.9999204 

1.00310 

443 

1.6133271 

1.6133275 

1.0000002 

0.991158 

490 

1.1106925 

1.1106931 

1.0000005 

0.959938 

510 

0.65016730 

0.65016802 

1.0000012 

0.985839 

555 

0.27024732 

0.27024714 

0.99999934 

0.993857 

670 

0.014412810 

0.014454512 

1.0069607 

0.959650 

765 




0.946 

865 




1.000 


of these other estimators had a negligible effect (<0.1%) 
on the vicarious gains. 

The vicarious calibration of band 6 is difficult because 
of the low values of Lwn in this band. Several of the Sea- 
WiFS scenes yield negative Lwn f° r band 6. Additionally, 
the calibration of MOBY over this band pass is problematic 
because of the crossover between the two spectrometers (D. 
Clark, pers. comm.). As a result, the matchups where neg- 
ative Lwn values occurred for band 6 were excluded from 
the vicarious calibration. 

Figure 31 shows a number of matchups that are out- 
liers. These outliers are defined to be individual match-up 
data points for any of bands 1-5 that deviate from the 
mean ratio in that band by more than two standard de- 


viations (2a). The CVT has attempted to determine why 
the outlier matchups plotted in these figures deviate so far 
from the mean L%$ ratios. The CVT has looked for corre- 
lations between the outliers and solar zenith angles, space- 
craft zenith angles, aerosol optical depth, e value, chloro- 
phyll concentration, local wind speed, and ozone concen- 
tration without success. The use of the geometric mean in 
estimating the L$$ ratios allows the outlier matchups to 
be included in the vicarious calibration. 

The match-up time series plotted in Fig. 31, excluding 
the outlier matchups, do not show any trends with time, 
which indicate the time corrections applied to bands 1,2, 
5, 6, 7, and 8 do not have any significant residual errors. 
Earlier during the mission, a 0.5% error in the ratio of 
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bands 7:8 of gave rise to a discernable error in the match- 
up time series. The matchups will be discussed as functions 
of scan angle in Eplee and McClain (2000). 

6.4 DISCUSSION 

The departure of the vicarious gains for bands 1-6 from 
unity are the result of uncertainties in the atmospheric 
correction algorithm, uncertainties in the laboratory cal- 
ibration of SeaWiFS, and uncertainties in the laboratory 
calibration of MOBY. The uncertainty in the SeaWiFS cal- 
ibration, as estimated by the vicarious gains, ranges from 
1-4% depending on the band. This result is consistent with 
the uncertainty derived from the prelaunch recalibration of 
SeaWiFS (Johnson et al. 1999) and with the uncertainty 


derived from the calibration transfer-to-orbit experiment 
(Barnes et al. 1999b). 

The time series of match-up values provides the 
CVT with a check on the time corrections applied to the 
individual bands. Currently, the time series show no signif- 
icant residual errors. As the number of matchups increases 
with time, the accuracy of these checks will increase ac- 
cordingly. 

The latest set of vicarious gains should be applied to 
the SeaWiFS data in conjunction with the current Sea- 
WiFS calibration table, which contains the temporal cor- 
rections for the instrument. The CVT will periodically up- 
date the vicarious calibration coefficients over the five-year 
SeaWiFS mission, e.g., at each reprocessing, as additional 
matchups become available. 
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Chapter 7 


The Updated SeaWiFS Coccolithophore Algorithm 

Wayne D. Robinson 
SAIC General Sciences Corporation 
Beltsville, Maryland 

Abstract 

A coccolithophore identification algorithm has been adapted to work reasonably well for SeaWiFS data contain- 
ing coccolithophore blooms at a number of sites in the Bering Sea and the North Atlantic. Changes were made in 
the initial algorithm thresholds that resulted in considerable improvements in the detection of coccolithophores 
for some cases. Future studies should include additional conditions which were considered in the formulation 
of the original coccolithophore algorithm, but not considered here — red tides, sediments, whitings, and haze. 
Finally, a future study may want to re-examine the predictors that are employed to detect the coccolithophore 
signal. The SeaWiFS 443, 510, and 555 nm bands are slightly different than the CZCS bands used to derive the 
current predictors. In addition, SeaWiFS has bands at 412 and 490 nm, which CZCS did not have, and these 
might be useful in deriving a better set of predictors to detect coccolithophores and other water types. 


7.1 INTRODUCTION 

The algorithm for detecting coccolithophores in Sea- 
WiFS data was initially developed with the aid of CZCS 
measurements taken for various water conditions, includ- 
ing the presence of coccolithophores (Brown and Yoder 
1994, and Brown 1995). The algorithm was developed by: 

a) Empirically determining the spectral signatures of 
coccolithophore blooms and various common non- 
bloom conditions; 

b) Deciding on feature predictors, i.e., spectral quan- 
tities which are correlated to the presence of coccol- 
ithophores, to be used in the algorithm; and 

c) Establishing decision boundary values for each pre- 
dictor which would allow the blooms to be spec- 
trally distinguished from other conditions. 

Predictors were derived from the normalized water- 
leaving [Lwn) radiances in three CZCS bands: 440, 550, 
and 520 nm and their ratios. Hereafter, the Lwn ratio 
constructed from band A divided by band B is indicated 
by L£, the Lwn(490)/Lww( 555) band ratio is rep- 
resented by L555. The following five predictors were em- 
ployed in the classification algorithm: 

1. Lwn (440), 

2. £wn(550), 

0 r 440 
^5201 

4. i« 8 , 

5. I'slo? an d 


6. L a (670) (Brown and Yoder 1994). 

These predictors were incorporated into the SeaWiFS 
level -2 processing program to derive the coccolithophore 
flag for each image pixel (McClain et al. 1995 and Darzi 
1998). When the first SeaWiFS-observed coccolithophore 
bloom was found in the Bering Sea (Vance et al. 1998), it 
was discovered that the algorithm was doing a poor job 
of identifying the bloom, which could clearly be seen in 
the true color imagery. The predictors were adjusted so 
that this case, and others during the same bloom, would 
be identified correctly. 

In the months of April-Juiy 1998, many coccolitho- 
phore bloom cases were observed in the North Atlantic 
and the Bering Sea. Comparisons of the coccolithophore 
flagged areas to the imagery showed that the algorithm 
was still not identifying part or all of these blooms. 

The availability of more coccolithophore test cases and 
the poor performance of the algorithm on some of these 
scenes has prompted the re-examination of the coccolitho- 
phore algorithm. The work done to improve the algorithm 
is described here. 

7.2 THE COCCOLITHOPHORE TEST 

The coccolithophore detection algorithm was developed 
to identify coccolithophore blooms and separate their sig- 
nal from those of clear water, red tide, haze, and sedi- 
ments. Although it was hoped that the algorithm could 
distinguish whitings — suspended lime muds — from cocco- 
lithophores, the previous work (Brown 1995) was unable 
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to distinguish them at all. In addition, the algorithm had 
difficulty with sediment and atmospheric haze. 

The coccolithophore identification tests shown in Ta- 
ble 12 use the aforementioned predictors to eliminate all 
conditions not caused by coccolithophores. Any remaining 
data are classified as coccolithophores. 

Table 12. The test sequences to identify coccol- 
ithophores. The center column indicates the test 
on the SeaWiFS parameter or ratio, and the right 
column indicates the condition that exists if this 
test is true. If none of tests 1-6 are true, the data 
is considered to contain a significant amount of coc- 
colithophores. 


No. 

Test 

Condition(s) 

1 

Lwn{ 550) < 0.8 

Clear water 

2 

Lvw(440) <1.1 

Clear water and 
possibly sediments 

3 

£ 55 ° < 1.0 

Sediments and red tide 


£55“ > 2.0 

Haze 

4 

^520 ^ 0.95 

Sediments 


^520 >1-5 

Haze 

5 

£555 < 1.0 or 

Unknown 


LUi > 16 

condition 

6 

L o (670) > 1.1 

Haze 

7 

Remaining points 

Coccolithophores 


Several comments should be made about the tests ap- 
plied to the CZCS data and the thresholds that were used. 
The test with L o (670), the aerosol radiance at 670 nm, was 
added by Brown and Yoder (1994) as a test for haze con- 
tamination and was important for use with CZCS data 
because the CZCS atmospheric correction assumed a fixed 
aerosol type over an entire scene. This test should now 
be less critical as SeaWiFS processing corrects for haze on 
a pixel-by-pixel basis. The bands used for SeaWiFS are 
slightly different — 443, 510, and 555 nm, with the corre- 
sponding normalized water-leaving radiances of Lvwv(443), 
Lwn{ 510), and Twjv(555) — as compared to 440, 520, and 
550 nm for CZCS (see Hooker et al. 1992 for CZCS and 
SeaWiFS descriptions). 

In addition to the central wave numbers being differ- 
ent, the response functions of the bands are not the same 
for CZCS and SeaWiFS. To detect coccolithophores with 
SeaWiFS, it could be expected that some changes would 
have to be made, at least in the test thresholds if not in the 
predictors. It should also be expected that in comparisons 
of coccolithophore-laden waters with clear water, the first 
two tests (Table 12), which are mainly concerned with dis- 
criminating coccolithophores from clear water, should be 
the determining tests. 

The initial tests of this algorithm with SeaWiFS data 
were done in the first month of operation. At that time, 
very little data were available containing coccolithophores. 


When the algorithm was applied to an observed coccolith- 
ophore bloom in the Bering Sea (Vance et al. 1998) from 
18 September to 3 October 1997, it was found that the 
above tests completely missed the identification of cocco- 
lithophores and identified neighboring clear water as hav- 
ing coccolithophores. It was assumed that the predictors 
remained valid and the test thresholds required adjust- 
ment. 

By examining points in the clear and coccolithophore 
areas, the tests that determine a coccolithophore bloom 
were set to the following: 

if Lvvw(555) > .81 and Z/Wjv( 443) >1.1 and 
L a (670) < 1.1 and 

0.60 < Liil < 0.90 and (24) 

0.90 < LlH < 1.20 and 
0.60 < L51Q < 0.90, 

then this is a coccolithophore pixel. 

These were relatively radical changes, but they iden- 
tified the coccolithophore bloom in the Bering Sea well 
(Figs. 32a— c), and served as the initial coccolithophore al- 
gorithm. With a lack of other cases at the time, these 
tests were adopted for operational processing as the best 
set possible. 

7.3 TESTS ON NEW CASES 

In the months of April- July 1998, a number of new 
coccolithophore cases appeared, especially in the North 
Atlantic. Poor performance of the mask in these Atlantic 
blooms and the availability of more test cases prompted 
these latest tests to improve the coccolithophore algorithm. 

In this study, three cases are examined to aid in the 
improvement of the coccolithophore algorithm. 

1. 20 September 1997 at 0014 UTC: This case is actu- 
ally one of the original cases in the Bering Sea that 
was used to tune the thresholds for SeaWiFS data; 
it will be useful here to verify that the changes still 
apply to the previous cases. 

2. 6 April 1998 at 1629 UTC: This is a coccolithophore 
bloom which occurred in the North Atlantic (57° W, 
40°N). 

3. 20 July 1998 at 0046 UTC: This is another bloom 
in the Bering Sea. 

Images are shown of the first two cases (Figs. 32 and 33) 
to show the performance of the initial and revised algo- 
rithms. 

Figure 32 shows the 20 September 1997 case. Fig- 
ure 32a shows an image of the 443 nm level-la data. This 
band was chosen because it does the best job of displaying 
the dark ocean, the bright clouds, and the hazy appear- 
ance of the coccolithophore signal. In true-color images, 
made by assigning the 670, 555, and 412 nm bands to red, 


52 




McClain, Ainsworth, Barnes, Eplee, Patt, Robinson, Wang, and Bailey 



Fig. 32 . Images of the coccolithophore bloom seen by the SeaWiFS 443 nm band in the Bering Sea on 
20 September 1997 at 0014 UTC. a) The 443 nm radiances and an overlay of the coastlines of Alaska 
and Siberia are shown in black, b) The same band is shown with the grey scale inverted to make the 
coccolithophore area appear dark (the masked areas of clouds, land, and high radiances areas are shaded 
white, leaving the coccolithophore signal as a dark area in the center), c) The same scene as in panel 
b is given, but with the coccolithophore bloom identified by the initial coccolithophore identification 
algorithm in black, d) The coccolithophore is identified in black using the revised algorithm. 
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Fig. 33. The same as Fig. 32, but for a region of the North Atlantic at (57°W,40°N) on 6 April 1998 
at 1629 UTC. 


green, and blue respectively, this appears as a distinctly 
aquamarine region. Figure 32b is an enhanced version of 
Fig. 32a showing the coccolithophore signal as a darker 
area in the center of the image. The areas masked for 
land, high radiances, and clouds are shaded white. The 
masked region in the center of the coccolithophore bloom 
(white region) was caused by the unusually high radiances 
of the bloom. The radiances were so high that they were in 
the poor resolution portion of the SeaWiFS instrument’s 
detector response (that is, above the knee in the instru- 
ment bilinear response), which is a masking condition. In 
Fig. 32c, the pixels the initial algorithm identified as hav- 
ing coccolithophores have been shaded black on a white 
background. It can be seen that the initial coccolithophore 
algorithm appears to handle most of the pixels containing 
coccolithophores (less the masked areas). This would be 
expected as this was one of the initial training sites. 

A similar treatment is shown for the coccolithophore 
bloom in the Atlantic Ocean that occurred on 6 April 1998 
(Figs. 33a-c). For this bloom, only a small fraction of the 
coccolithophore- affected pixels are identified by the initial 
algorithm (Fig. 33c). The coccolithophore flag did not 
work well on the 20 July 1998 case either, although it per- 
formed better than in the 6 April case. 


Each of the six predictors were examined using selected 
pixels from the three cases above to determine why the coc- 
colithophore test was not working well on the other cases. 
The six predictors examined were Lwn{ 443), Lww(555), 
L o (670), L555, and For each case, groups of 

points were selected from three regions of the images with 
about three observations per region. The first region con- 
tained pixels that were not affected by coccolithophores 
in clear water. The second region contained pixels that 
had visible coccolithophores, but were not identified by 
the algorithm as having coccolithophores, and the third 
had pixels correctly identified as having coccolithophores. 

Figures 34a-f display each of the predictor values for 
the points from the three regions for the three cases con- 
sidered. To aid in identification, the points are spread 
in the horizontal and marked with different symbols to 
separate the three regions: clear (with a diamond sym- 
bol); coccolithophore, but not identified (asterisk symbol); 
and coccolithophore identified (box symbol). The points 
are also spread horizontally to segregate observations from 
each of the three cases. In addition to the pixel values, solid 
horizontal lines are plotted to show the initial thresholds 
determined for SeaWiFS. 
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Fig. 34. Graphical plots of the six coccolithophore tests and SeaWiFS test pixels taken in clear waters 
(diamonds), coccolithophore waters that were not identified by the initial algorithm (asterisks) and 
coccolithophore waters that were identified by the initial algorithm (squares). The points were taken 
from three separate coccolithophore events: 20 September 1997 at 0014 UTC in the Bering Sea, 6 April 
1998 at 1629 UTC in the North Atlantic, and 20 July 1998 at 0046 UTC in the Bering Sea. The solid 
horizontal lines indicate the initial test thresholds, while the dashed lines in panels d, e, and f show 
the revised high test thresholds. Each test for coccolithophores is positive in the range spanned by the 
arrow on the right side of the plot. The tests are for: a) Lww(443), b) Lww( 555 )< c) L„(670), d) L 555 , 

e ) -^555’ an< ^ 0 ^5 10' 
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The tests for L WN { 443) and L WN { 555) (Figs. 34a-b) 
are functioning properly for all of the coccolithophore iden- 
tified cases and for two of the three coccolithophore not- 
identified cases. The exception is the 20 September 1997 
case, for which the coccolithophore identification already 
works well. The misidentified points from this case will 
not be considered as strongly in making revised thresholds 
because they could represent more complex conditions. 

The test for L a (670) (Fig. 34c) also appears to be func- 
tioning well although it is clear that it is of little use in 
discriminating most of the pixels chosen. This is under- 
standable considering the L a (670) test is used to remove 
pixels with high aerosol; the tests in this study use either 
clear or coccolithophore containing pixels. 

The T555, T555, and tests (Figs. 34d-f) failed to 
identify many of the misidentified coccolithophore pixels. 
In each of these tests, the misidentified pixels have val- 
ues greater than the high threshold which w as initially set. 
As with the L a (670) test, these tests show little discrimi- 
nation between the dear-water and coccolithophore pixels 
because these tests are used to identify conditions that are 
not included in the present study (i.e., haze, red tides, and 
sediments) . 

High thresholds in the tests for Ljg, and L\\l 

are adjusted by consulting these plots, applying the values 
to the test cases, and then viewing the resulting masked 
pixels. The dotted horizontal line in Figs. 34d-f represent 
the value of the revised high thresholds. 

The coccolithophore test with the modified thresholds 
becomes: 

if Lvwv(555) > .81 and Lw7sr(443) > 1.1 and 
L a (670) <1.1 and 

0-60 < L555 <1.10 and (25) 

0.90 < Lf*° < 1.32 and 

0.60 < < 0.92, 

then this is a coccolithophore bloom pixel. 

The thresholds were increased in a conservative manner 
for the Lggg, L555, and tests to correctly identify pix- 
els in the test cases. These thresholds could be expanded 
more, but that should be done only when haze, sediment, 
and red tide test cases are included in the study. 

Figures 32d and 33d show the results of the revised 
thresholds which do little to change the distribution of 
identified coccolithophores in the 20 September 1997 case 
(Fig. 32d), because the initial thresholds were properly 


identifying most of the coccolithophore pixels. The 6 April 
1998 case (Fig. 33d) however, showed a significant improve- 
ment in coccolithophore identification. The revised thresh- 
old also improved the 20 July 1998 case. The present algo- 
rithm does not appear to misidentify clear water as having 
coccolithophores. 

7.4 CONCLUSIONS 

The coccolithophore identification algorithm of Brown 
(1995) has been adapted to work reasonably well for Sea- 
WiFS data containing coccolithophore blooms at a number 
of sites in the Bering Sea and the North Atlantic. In some 
of the newer cases found, the improvements in the test- 
ing thresholds have significantly increased the ability to 
detect coccolithophore blooms, although many aspects of 
the testing can be improved. 

The coccolithophore algorithm was tuned to one case 
in the Bering Sea and subsequently improved with new 
cases; however, the number of cases needs to be increased 
in number and geographical location. This can be accom- 
plished by searching the present data for more coccolitho- 
phore cases and by getting more cases from the new data. 
This study used only a limited number of observations to 
adjust the thresholds. Future studies may benefit from the 
use of a greater number of coccolithophore and clear- water 
observations. 

The tests in this study only considered coccolithophore 
and clear-water conditions. This means that four out of 
the six testing conditions of the coccolithophore algorithm 
have not been investigated fully. Future studies should in- 
clude additional conditions which were considered in 
Brown’s study — red tides, sediments, whitings, and haze. 
The SeaWiFS Project is starting to compile examples of 
some of these other conditions.! 

Finally, a future study may want to re-examine the pre- 
dictors that are employed to detect the coccolithophore sig- 
nal. The SeaWiFS 443, 510, and 555 nm bands are slightly 
different from the CZCS bands used to derive the current 
predictors. In addition, SeaWiFS has bands at 412 and 
490 nm, which CZCS did not have, and these might be 
useful in deriving a better set of predictors to detect coc- 
colithophores and other water types. 
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f Browse the calibration and validation image gallery web site 
at: http : / / calval-2 . gsf c . nasa . gov/calval/. 
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Abstract 

Modifications to the atmospheric correction algorithm for the third SeaWiFS reprocessing are described in this 
chapter. The updates include changes to the aerosol look-up tables, the atmospheric diffuse transmittance 
tables, the ocean whitecap computations, and the implementation of new Rayleigh radiance tables generated 
with a variety of ocean surface wind speeds. In addition, computation of a new SeaWiFS atmospheric product, 
the Angstrom exponent, is described. These modifications significantly improve the SeaWiFS retrieval results. 


8.1 INTRODUCTION 

It is well known that atmospheric correction removes 
more than 90% of sensor-measured signals contributed from 
the atmosphere in the visible spectrum, and is the key pro- 
cedure in ocean color imagery data processing. With the 
successful launch of SeaWiFS (McClain et al. 1998) on 1 
August 1997 and its data processing since then, it is very 
important to periodically re-evaluate and update the Sea- 
WiFS atmospheric corrections (Hooker et al. 1992). 

The SeaWiFS atmospheric correction algorithm uses 
two near-infrared (NIR) bands (765 and 865 nm) to esti- 
mate the aerosol optical properties which are then extrap- 
olated into the visible part of the spectrum (Gordon and 
Wang 1994a). The implementation of the algorithm was 
achieved using look-up tables for Rayleigh scattering, aero- 
sol contributions, and the effects of the atmospheric diffuse 
transmittance. Ocean whitecap contributions at the Sea- 
WiFS bands were estimated using a reflectance model. In 
this chapter, a brief overview of the SeaWiFS atmospheric 
correction algorithm is presented, and outlines are given 
of some updates and modifications in the aerosol look-up 
tables, the atmospheric diffuse transmittance tables, and 
computations of the whitecap contributions. 

8.2 ALGORITHM 

To better describe the SeaWiFS atmospheric correction 
algorithm and its implementation into the data process- 
ing system, the reflectance p = ttL/(Fq cos#o) is defined, 
where L is the radiance in a given solar and viewing geom- 
etry, Fq is the extraterrestrial solar irradiance, and 6o is 
the solar zenith angle. The SeaWiFS measured reflectance 


at the TOA can be written as: 

Pt(X) — Pr(X) + Aj(A) + Pra{X) (26) 

+ t(X)pf( X) -f f(A)pvv(A), 

where p r (A), p a (A), and p ro (A) are the contributions from 
multiple scattering of air molecules (Rayleigh scattering), 
aerosols, and Rayleigh-aerosol interactions, respectively; 
p/(X) is the reflectance at the sea surface which arises from 
sunlight and skylight reflecting from foam (or whitecaps) 
on the surface (Gordon and Wang 1994b); and pw{X) is 
the water-leaving reflectance, which is the desired quantity 
in ocean color remote sensing. The t(X) term is the atmo- 
spheric diffuse transmittance (Wang 1999, and Yang and 
Gordon 1997) which accounts for the effects of propagat- 
ing water-leaving and whitecap reflectances from the sea 
surface to the TOA. The surface sun glint term has been 
ignored in (26). 

Because more than 90% of the signal in the visible spec- 
trum measured at the satellite altitude is contributed by 
the atmosphere and ocean surface effects — the first four 
terms in (26) — accurately removing these effects is crucial 
to the success of any ocean color remote sensing exper- 
iment. The SeaWiFS atmospheric correction algorithm 
(Gordon and Wang 1994a) uses the two SeaWiFS NIR 
bands centered at 765 and 865 nm to estimate the atmo- 
spheric effects and extrapolate these into the visible spec- 
trum. Unlike Rayleigh scattering, which can be computed 
accurately, aerosol scattering is highly variable, and the 
effects of pa{ X) + p ra (A) in (26) on the imagery cannot 
be predicted a priori . The water-leaving reflectance at the 
two NIR bands, however, is usually negligible because of 
strong water absorption. The radiances measured at these 
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two NIR bands, therefore, are essentially the contributions 
from the atmosphere. For the two SeaWiFS NIR channels, 
(26) can be written as 

p t (X) - Pr(A) - t(\)p f {\) = p a { A) 4- p ro (A). (27) 

The effects of p a (A) *f p ra (A) interactions in the imagery, 
therefore, can be estimated at the two NIR bands from the 
sensor-measured radiances, the computed Rayleigh scat- 
tering reflectances, and the estimated whitecap contribu- 
tions (Gordon and Wang 1994b). This quantity is then 
extrapolated and removed in the visible. 

The extrapolation was achieved through a process of 
aerosol model selection based on an evaluation of the at- 
mospheric correction parameters, e(A*,Aj). The latter is 
defined as (Gordon and Wang 1994a, and Wang and Gor- 
don 1994): 

e(A ^> ■ £$■ < 28 > 
where p 0 «(Ai) and Paa(Aj) are the single scattering aerosol 
reflectances at wavelengths A* and A ,, respectively. The 
A j value is usually taken as the longest NIR band (865 nm 
for SeaWiFS). The e(A*,Aj) value characterizes the spec- 
tral variation of the aerosol extinction coefficient which 
includes the aerosol optical thickness, single scattering al- 
bedo, and aerosol scattering phase function, thereby link- 
ing the value of e(A*, \j) and the aerosol model. 

The implementation of the atmospheric correction al- 
gorithm into the SeaWiFS data processing system was 
achieved through the use of look-up tables based on a 
large number of radiative transfer simulations (approxi- 
mately 25,000) which use the aerosol models developed by 
Shettle and Fenn (1979). The main look-up tables con- 
tain information on the p a (A) + Pro (A) values for various 
aerosol optical and microphysical properties (for a number 
of different aerosol models with a variety of aerosol opti- 
cal thicknesses) along with solar and viewing geometries 
at the eight SeaWiFS spectral bands. The two other ta- 
bles, which are much smaller in size, are the p r (A) and *(A) 
tables. 

To obtain the value of pw{ A) in (26), four quanti- 
ties have to be estimated: a) p r (A), b) the reflectance of 
Pa (A) + Pra(A), c) £(A), and d) p/(A). In the following four 
sections, some updates are discussed for the aerosol look- 
up tables Po(A) + p ra (A), the atmospheric diffuse transmit- 
tance tables, modifications in computing e(A t , A j) when its 
value lies outside of the model range, and modifications in 
computing the foam reflectance contribution. 

8.3 AEROSOL LOOK-UP TABLES 

Before the reprocessing (prior to August 1998), the Sea- 
WiFS aerosol look-up tables were generated using 12 aero- 
sol models .with 8 aerosol optical thicknesses (0.05, 0.1, 
0.15, 0.2, 0.3, 0.4, 0.6, and 0.8) and various solar and view- 
ing geometries. The 12 aerosol models are based on three 


aerosol types, maritime, coastal, and tropospheric, and 
four relative humidity (RH) values of 50, 70, 90, and 99. 
The maritime and tropospheric aerosol models are from 
Shettle and Fenn (1979), whereas the coastal model is de- 
rived from their works as described in Gordon and Wang 
(1994b). These aerosol models are believed to be repre- 
sentative of the aerosol optical properties over the ocean. 
Table 13 provides the 12 aerosol model names, model num- 
bers, and their corresponding RH values used before the 
reprocessing. These 12 aerosol models were used as can- 
didates for generating the aerosol look-up tables for the 
atmospheric corrections. 

Table 13. The 12 aerosol models used before the 
reprocessing. The model names are formed from a 
one-letter code indicating the model type, “C” for 
coastal, “M” for maritime, and “T” for troposheric, 
followed by a two-digit RH value. 


No. 

RH [%] 

Names 

1-4 

50, 70, 90, and 99 

M50 . 

. M99 

5-8 

50, 70, 90, and 99 

C50 . 

. C99 

9-12 

50, 70, 90, and 99 

T50 . 

. T99 


Figures 35a and 35b provide examples of e(A,865) as 
a function of wavelength for these 12 aerosol models. Fig- 
ure 35a is for the case of 0q = 0° and a sensor zenith 
angle of 0 = 45°, whereas Fig. 35b is for 0 0 = 60° and 
0 - 45 0 with a relative azimuthal angle of <f> = 90°. As 
discussed in the previous section, the e(A,865) values, as 
shown in Figs. 35a and 35b, were used as guides to select 
the aerosol model, and its optical and radiative properties 
were extrapolated from the SeaWiFS NIR bands into the 
visible spectrum in the atmospheric corrections. 

As shown in Figs. 35a and 35b, the T50 and M99 aero- 
sol models give the highest and lowest e(A,865) values 
within these 12 aerosol models, respectively. The T50 and 
M99 aerosol models, therefore, were taken as default mod- 
els such that, if the retrieved e(765,865) value is higher 
(lower) than the value predicted by the T50 (M99) model, 
the retrieval model is defaulted to the T50 (M99) aerosols. 
Figures 35a and 35b show that in the single scattering 
case, for the M99 model, the aerosol reflectance at 412 nm 
contributes a factor of approximately 0.85-0.93 at 865 nm, 
while for the T50 model, the aerosol reflectance at the 
412 nm wavelength is a factor of approximately 2. 2-2. 4 
higher than at the 865 nm wavelength. 

In the SeaWiFS data processing, however, there was 
evidence that some new aerosol models needed to be in- 
cluded, in particular, for the optical properties with lower 
e(A, 865) values: 

1. It was found that there was a significant number 
of pixels defaulting to the M99 aerosol models, in- 
dicating that the lowest e(765, 865) value predicted 
by the M99 model was still too high. For example, 
some case studies showed that about 35-40% of the 
retrievals defaulted to the M99 aerosol model. 
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Fig. 35. The c( A, 865) value as a function of wavelength for the 12 aerosol models used before the reprocessing 
for the sensor zenith angle of 6 = 45°, and the solar zenith angle of a) 0 O = 0° and b) 0 O = 60°. 
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2. It appeared that SeaWiFS retrieved more pixels as- 
sociated with the negative water-leaving reflectance 
at the short wavelengths than at the long wave- 
lengths, which indicated that probably for some 
cases, the aerosol contributions at the short wave- 
lengths were overestimated. 

For these reasons, the decision was made to update the 
aerosol look-up tables. 

Following a reccommendation from H. Gordon (pers. 
comm.), the oceanic aerosol model from Shettle and Fenn 
(1979) was added into the candidate models. During the 
second reprocessing (August 1998 to December 1999), the 
C70 and T70 models were replaced by the oceanic aerosols 
with an RH of 90 and 99% (090 and 099, respectively). 
It was found, however, that in some cases (depending on 
the solar and viewing geometry), the e values predicted by 
the 090 model intersected with that of the M99 model (de- 
pending on the wavelength), e.g., the 090 model predicts a 
higher e value at 765 nm and a lower value at 443 nm than 
those of the M99 model which leads to discontinuities in 
the SeaWiFS derived products. 

In the third reprocessing, therefore, only the T70 model 
was replaced by the 099 model. The number of aerosol 
models was kept at 12, and Table 14 shows the 12 aerosol 
models used in the updated aerosol look-up tables. Us- 
ing the same solar and viewing geometries as in Fig. 35, 
Figs. 36a and 36b provide the values of e(A, 865) as a func- 
tion of wavelength, A, for the updated aerosol models. 


Table 14. The updated 12 aerosol models used 
in the reprocessing which include the new oceanic 
model type (099). 


No. 

RH {%} 

Names 

1 

99 

099 


2-5 

50, 70, 90, and 99 

M50 . 

. M99 

6-9 

50, 70, 90, and 99 

C50 . 

. C99 

10-12 

50, 90, and 99 

T50 . . 

. T99 


With the Oceanic aerosol model, the lowest e(A,865) 
value is expanded. The 099 model predicts a factor of ap- 
proximately 0.7 single scattering aerosol reflectance con- 
tribution at 412 nm as at 865 nm; this value is lower than 
the value from the M99 model (approximately 0.85-0.93). 
Inclusion of the new models decreased the number of pix- 
els being processed with the default aerosol model by more 
than one half. 

In addition, the computation of the e(A, 865) values was 
modified for cases where the retrieved e(765, 865) value was 
lower than the lowest predicted values of the 12 models 
(i.e., the 099 model). Instead of using the default model to 
compute e(A,865), an analytical formula (Wang and Gor- 
don 1994) was used with retrieved e(765, 865) values to 
estimate e(A, 865): 


e(A, 865) = exp 


In (6(765,865)) 


865 - A 
100 


(29) 


With these changes, the retrieved results are signifi- 
cantly improved by using the updated aerosol look-up ta- 
bles (Robinson et al. 2000). 

8.4 TRANSMITTANCE TABLES 

Because the t(A) values depend on the aerosol mod- 
els, it is necessary to update the transmittance tables with 
the updated aerosol models, i.e., including the transmit- 
tance tables for the 099 model. In the old tables (prior to 
August 1998), there was an error in the atmospheric dif- 
fuse transmittance tables, in which the Fresnel- reflecting 
ocean surface was mistakenly not included in the compu- 
tations. This caused the SeaWiFS-retrieved normalized 
water-leaving radiance to depend strongly on the SeaWiFS 
scan angle, in particular, at the SeaWiFS scan edges. This 
error was corrected in the updated diffuse transmittance 
tables, and the reprocessing results show the normalized 
water-leaving radiances now have no obvious dependence 
on the SeaWiFS scanning angles (R. Eplee, pers. comm.). 

8.5 WHITECAP CONTRIBUTIONS 

As discussed in Sect. 8.2, the foam contribution in the 
SeaWiFS imagery is estimated by using a reflectance model 
(Koepke 1984) based on the sea-surface wind speed (Gor- 
don and Wang 1994b). It is assumed that the foam is 
white, i.e., the whitecap reflectance is independent of wave- 
length. The in situ measurements, however, show a signif- 
icant uncertainty in p/(X) with the sea-surface wind speed 
(e.g., Monahan 1971). The measurement data are par- 
ticularly noisy for sea-surface wind speeds greater than 
7-8 ms* 1 wherein the uncertainty is usually greater than 
100% (Fig. 1 in Gordon and Wang 1994b). 

In recent studies, both Frouin et al. (1996) and Moore 
et al. (1998) found that contrary to the previous measure- 
ments, the whitecap reflectance is spectrally dependent. 
The reflectance contributions are significantly smaller at 
the near-infrared than in the visible spectrum, because of 
the stronger ocean water absorption at the longer wave- 
lengths. This reduces the reflected photons from the sub- 
merged bubbles. The SeaWiFS observations also show that 
the whitecap reflectance model, used in the atmospheric 
corrections, likely overestimated the whitecap contribu- 
tions significantly, in particular, for the larger sea-surface 
wind speed. The overcorrection was particularly evident 
in very clear open ocean scenes with high surface wind, 
which resulted in low normalized water-leaving radiance 
retrievals. 

The SeaWiFS Project has adopted the Frouin et al. 
(1996) results (with updates) for the spectral dependence, 
and the results from Moore et al. 2000) for the magnitude 
of the whitecap contributions. This leads to a reduction 
in the whitecap radiance contributions for SeaWiFS bands 
1-8 by factors of 0.4, 0.4, 0.4, 0.4, 0.4, 0.3557, 0.3040, and 
0.2580, respectively, relative to the values estimated us- 
ing Gordon and Wang (1994b). The whitecap radiance for 
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Fig. 37. The e(A,865) value as a function of wavelength A, for the SeaWiFS updated 12 aerosol models. 


wind speeds greater than 8 ms" 1 are set equal to the value 
computed at 8ms -1 . Results from in situ measurements 
showed similar characteristics (Moore et al. 2000). It is 
believed, however, that some validation efforts with more 
experimental data are needed to further fine tune the foam 
reflectance model. In addition, evaluations of the National 
Center for Environmental Prediction (NCEP) wind speeds 
using buoy winds near Hawaii and Bermuda show very 
good agreement, and differences between SeaWiFS and the 
Bermuda Test Bed Mooring (BTBM) water-leaving radi- 
ances were uncorrelated to wind speed (S. Bailey, pers. 
comm.) which suggests the 8 ms -1 limit is not introduc- 
ing a bias in the SeaWiFS products. In fact, it is better 
to underestimate foam radiance than to overestimate it 
(Gordon and Wang 1994b), because the residual will be 
included in the aerosol radiance. 

8.6 NEW RAYLEIGH TABLES 

New Rayleigh radiance look-up tables for the SeaWiFS 
bands were generated using the method developed by Gor- 
don and Wang (1992) for the various ocean surface wind 
speeds. A bidirectionally shadowing factor for a collection 
of individual wind-ruffled facets (Gordon and Wang 1992) 
was used in all computations. These new wind speed de- 
pendent Rayleigh tables were implemented in the third re- 
processing. Prior to the third reprocessing, Rayleigh tables 
were generated with a flat ocean surface assumption (where 
the wind speed was zero). The Rayleigh radiance tables 
were generated at eight ocean surface wind speeds corre- 
sponding to 0, 1.9, 4.2, 7.5, 11.7, 16.9, 22.9, and 30.0 ms -1 , 


respectively. Tests showed that, with the new Rayleigh 
tables, the SeaWiFS-derived ocean color products are sig- 
nificantly improved, in particular, for cases of the larger 
solar zenith angles, i.e., 9q > 60°. A detailed study for the 
effects of the surface roughness on the SeaWiFS derived 
ocean color products is currently underway. 

8.7 Angstrom exponent 

The aerosol Angstrom exponent is widely used in the 
atmospheric aerosol and radiation communities. For two 
wavelengths at A* and A^, the Angstrom exponent is de- 
fined as 


Ta(Ai) 

r a(A j) 

A 1 q (A».A>) 

Ad 

(30) 

Qr(Aj,A^) = 

ln [sS$] 

l»fel ’ 

(31) 


where r a (Ai) and r a (Xj) are aerosol optical thicknesses 
measured at A, and Aj, respectively. The Angstrom ex- 
ponent, which is independent of the solar and viewing ge- 
ometry, can be used to relate the aerosol microphysical 
properties (particle size) and its optical spectral depen- 
dence. Figure 37 provides values of a(A, 865) as a function 
of wavelength, A, for the 12 SeaWiFS updated aerosol mod- 
els. Apparently, the a(A, 865) value is nearly independent 
of the wavelength, in particular for the oceanic, maritime, 
and coastal aerosol models. 
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With the two retrieved aerosol models in the atmo- 
spheric correction process, the SeaWiFS a(A T ,865) can be 
easily derived: 

a(Ai,865) = (1 — r) a a (A t , 865) + ra b (A t1 865), (32) 

where a°(A t ,865) and a b (A i5 865) are, respectively, the 
Angstrom exponent from derived models a and 6, and r 
is the ratio between two aerosol models from the retrieved 
c values. The SeaWiFS a(A t ,865) retrieval routine can be 


very easily implemented in the SeaWiFS data processing. 

It was decided that the a(510,865) data will be rou- 
tinely retrieved and archived. The reasons for choosing 
the 510 nm channel are that it is closely matched with the 
ground in situ measurement data (e.g., AERONET data) 
and it gives a good approximation to be used in deriving 
aerosol optical thickness at other SeaWiFS wavelengths 
(e.g., 412 nm). Nevertheless, the SeaDAS code can be eas- 
ily modified to obtain a(A t ,865) in any of the SeaWiFS 
wavelengths. 
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Chapter 9 


Correction of the Sun Glint Contamination on the SeaWiFS 
Aerosol Optical Thickness Retrievals 

Menghua Wang 

University of Maryland Baltimore County 
Baltimore, Maryland 

Sean W. Bailey 
Futuretech Corporation 
Greenbelt, Maryland 

Abstract 

For ocean color remote sensing, the measurement of radiances affected by sun glint has to be avoided or masked 
out. SeaWiFS has a capability of operationally tilting the sensor 20° away from nadir to minimize sun glint 
contamination, however, sun glint is still a factor near the subsolar point. In this chapter, results are presented 
which quantify the effect of sun glint contamination on the retrievals of atmospheric and bio-optical oceanic 
products. It was found that, although the sun glint contamination has a minor effect on the retrieved bio-optical 
oceanic products, the effect on the retrieved atmospheric products (e.g., aerosol optical thickness) is significant. 
A sun glint correction scheme is described, which was implemented in the SeaWiFS data processing. It was 

found that the sun glint correction significantly improves the derived atmospheric products in the vicinity of 
the subsolar point. 


9.1 INTRODUCTION 

In ocean color remote sensing, the radiance measured 
at the top of the ocean-atmosphere system can be written 
as stated in (13), 

L t {\) = L r { A) + L a { A) 

+ L ra ( A) + T(X)L g (X) (33) 

+ t(X)L f (X) + t(X)L w {X), 

where L r (A), L a ( A), and L ra (X) are the radiance contribu- 
tions from multiple scattering of air molecules (Rayleigh 
scattering), aerosols, and Rayleigh- aerosol interactions, re- 
spectively (Gordon and Wang 1994a). L g ( A) is the spec- 
ular reflection from the direct sun (sun glint) radiance, 
L/(X) is the radiance at the sea surface resulting from 
sunlight and skylight reflecting off foam (or whitecaps) on 
the surface (Gordon and Wang 1994b), and L w { A) is the 
water-leaving radiance. T( A) and t( A) are the atmospheric 
direct and diffuse transmittances at the sensor viewing di- 
rection, respectively. 

As there are usually no meaningful retrievals in the 
regions contaminated by sun glint, measurement of radi- 
ances affected by sun glint have to be avoided or masked 
out. SeaWiFS has the capability of operationally tilting 


the sensor 20° away from nadir to minimize sun glint con- 
tamination. Sun glint is still a factor, however, in the vicin- 
ity of the subsolar point. The SeaWiFS processing system 
computes L g ( A) from the Cox and Munk (1954) model 
based on the sea surface wind speed. A mask is applied 
to areas where the glint radiance is greater than a prede- 
termined threshold. Although the regions with the most 
significant sun glint contamination are masked out, there is 
still some residual sun glint contamination surrounding the 
mask. The sun glint contamination is particularly evident 
in the SeaWiFS derived atmospheric products, e.g., aerosol 
optical thickness. It is important, therefore, to develop a 
correction scheme for removing the effects of sun glint con- 
tamination. The ocean pigment concentration, however, is 
usually less affected by sun glint contamination, because 
the bio-optical algorithm uses a two-band ratio value in 
the derived water-leaving radiances (Fraser et al. 1997). 

In this chapter, the SeaWiFS sun glint mask from Cox 
and Munk (1954) is briefly described. Then, a sun glint 
contamination correction scheme and its implementation 
in the SeaWiFS data processing are proposed. Next, eval- 
uations of the NCEP wind speed data are presented, which 
are used in the SeaWiFS data processing in comparison 
with in situ measurements. Finally, some comparison re- 
sults are presented with and without sun glint corrections. 
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9.2 THE SeaWiFS SUN GLINT MASK 

It is convenient to rewrite the sun glint radiance as: 

Lg( A) = F 0 (A) To(A) Lgn, (34) 

where F 0 ( A) is the solar irradiance (adjusted for Earth-sun 
distance variations), T 0 ( A) is the atmospheric direct trans- 
mittance in the solar direction, and Lgn ls the normalized 
sun glint radiance, i.e. , the value of L g if there were no 
atmosphere and F 0 ( A) = 1. Note that for a given pixel, 
the Lgn value depends on the solar and viewing geometry, 
the sea surface wind speed, and the wind direction. 

The Lgn term is computed based on the assumption 
that the wind-roughened sea surface consists of a collection 
of individual facets obeying the Cox and Munk (1954) slope 
statistics. As an approximation, Lgn is computed with a 
further assumption that the wind-roughened surface slope 
distribution is independent of the wind direction. The in- 
puts for the computation are, at a pixel- by-pixel level, the 
solar and sensor viewing geometry, as well as the surface 
wind speed. The computed Lgn value is used as the sun 
glint mask and is applied to areas where the glint radiance 
is greater than a predetermined threshold. 

In the SeaWiFS data processing, the pixels for which 
Lgn > 0.005 are masked, and no further calculations for 
retrieving oceanic or atmospheric products occurs (Mc- 
Clain et al. 1995). There are also no corrections applied 
for pixels for which 0 < Lgn < 0.005, so it is likely 
that there will be some sun glint contamination outside 
the sun glint mask. It has been found that there is no 
apparent bias in the SeaWiFS retrieved oceanic products 
(e.g., water- leaving radiances and chlorophyll concentra- 
tion) around the subsolar point, but the derived aerosol 
optical thicknesses are always biased high in regions where 
0 < L G n < 0.005. To improve the SeaWiFS atmospheric 
products, sun glint contamination needs to be corrected. 

9.3 SUN GLINT CONTAMINATION 


The r r (A) and r 0 ( A) terms in (36) are the optical thickness 
for Rayleigh (air molecules) and aerosols, respectively, and 
Qo and 0 correspond to the solar and sensor zenith an- 
gles, respectively. It is assumed that the ozone absorp- 
tion effects in (35) and (36) have been removed. The sun 
glint corrected radiance, Lt, at the eight SeaWiFS bands 
can then be inserted into the atmospheric correction algo- 
rithm, replacing L t ( A), and data processing can proceed 
with corrected radiances. 

At the time of data processing, the aerosol optical thick- 
ness in (36) is unknown, so an iterative scheme is proposed. 
The sun glint corrected radiance, Lt, is first calculated with 
the uncorrected L*, wind speed, and an initial guess for r a 
[used here as r'(A0]. The initial L t is then used to derive 
r a which is then used as the input to derive the final sun 
glint corrected radiance and the final sun glint corrected 
aerosol optical thickness. 

Based on results from some case studies, the <(A t ) 
values are estimated using the following steps: 

a) Lgn is computed using the Cox and Munk (1954) 
model; 

b) Setting r a (865) = 0.1 (which is the global average 
value from SeaWiFS) in (36), La = L a + L ra at 
the SeaWiFS 865 nm band in (35) can be derived 
and then converted to the aerosol reflectance using 
pA = kLa/Fq cos a- n d 

c) A high <(A*) value is used for a low p A value (i.e., 
a <(865) with values of 1.0, 0.4, 0.2, and 0.12 are 
used, corresponding to p A of 0.001, 0.005, 0.008, 
and 0.01, respectively). 

The reason for using the large <(865) values for the low p A 
cases is to avoid the overcorrection of the sun glint contam- 
ination due to uncertainties in the Lgn estimation. These 
uncertainties arise from either uncertainty in surface wind 
speed or limitations of the Cox and Munk (1954) model, 
i.e., for the cases of small p A values (less than approxi- 
mately 0.005 corresponding to r a <« 0.05), confidence in 
the sun glint correction is low. 


If L g is known (34), it can be subtracted from the 
sensor-measured radiances (33), and subsequent data pro- 
cessing is based on the corrected radiances: 

L t = L r ( A) + L a (A) + L ra { A) 

4- t(X) Lf(X) + t(X)L w (X) y K 

where L t is the sun glint corrected radiance, which is de- 
fined as L t ( A) — Fo(A) Tq(A) T(A) Lgn > an< i 


T 0 (A)T(A) = exp 


-(r r (A) + r a (A)) 

/ 1 | 1 V 

\ cos 00 cos 0 / 


(36) 


9.4 WIND SPEED DATA EVALUATION 

To compute Lgn, the sea surface wind speed is needed. 
Based on a study and recommendation by Firestone et al. 
(1994), SeaWiFS uses the surface wind speed data pro- 
vided by NCEP which is gridded at 1°. The SeaWiFS 
data processing interpolates this coarse wind data to the 
SeaWiFS 1km pixel. Because the Cox and Munk (1954) 
model is sensitive to wind speed, some evaluations of the 
NCEP data were conducted. NCEP winds were compared 
to winds measured at the BTBM. Figure 38 provides an 
example scatter plot for the NCEP winds versus the in 
situ buoy data. The two data sets agreed reasonably well 
and there was no obvious bias of the NCEP wind data. 
This result lends confidence to the use of the interpolated 
NCEP winds in the Lgn calculation. 


65 









McClain, Ainsworth, Barnes, Eplee, Patt, Robinson, Wang, and Bailey 




Fig. 40. A quantitative comparison of the SeaWiFS derived r 0 (865) values for cases with and without sun 
glint contamination corrections: a) histograms from the results in Fig. 39, and b) data from pixels 70-460 at 
scan line 100 (dashed line in Fig. 39). 
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9.5 RESULTS 

To assess the efficacy of the proposed sun glint cor- 
rection scheme, some case studies were conducted. Fig- 
ure 39 provides an example of the SeaWiFS r„(865) im- 
ages for cases with and without sun glint corrections. Sea- 
WiFS data (file name S1998317034114) was acquired on 
13 November 1998 along the west coast of Australia at the 
location around -24° latitude and 120° longitude. The im- 
age in the top panel was processed without the sun glint 
correction, while the bottom image was processed with the 
correction applied. The value of r a (865) is scaled from 0- 
0.3. The sun glint mask (right side part of image, gray- 
white color) in Fig. 39 is clearly seen. In comparing two 
t 0 (865) images, it can be seen that the derived aerosol op- 
tical thickness is reduced around the area outside the glint 
mask when the correction is applied. 

Figures 40a and 40b show results of a quantitative 


comparison of these two images. Figure 40a compares the 
histograms from the two cases, while Fig. 40b compares a 
specific scan line (the dashed line in each image in Fig. 39). 
Figure 40a shows that, with the sun glint corrections, the 
large r 0 (865) values, which correspond to the sun glint 
contamination regions, are much reduced. The sun glint 
contamination effects sure clearly seen in Fig. 40b where the 
r 0 (865) is obviously biased, increasing as the pixel numbers 
increase (close to the glint mask region) for cases without 
the sun glint corrections. The T a (865) values are much 
more reasonable with the sun glint contamination correc- 
tion applied. 

The effects of sun glint contamination on the oceanic 
optical products were also studied. It was found that there 
is almost no effect on the SeaWiFS retrieved ocean prod- 
ucts, e g., normalized water-leaving radiances and chloro- 
phyll concentration. 
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Chapter 10 


Modifications to the TOMS Ozone 
Ancillary Data Interpolation 

Ewa J. Ainsworth and Frederick S. Patt 
SA 1C General Sciences Corporation 
Beltsville, Maryland 

Abstract 


The method for interpolating the TOMS ozone data used for level -2 processing was changed as a result of new 
information about the generation of the files by the TOMS Project. The new method and comparison of the 
old and new results are described in this report. 


10.1 INTRODUCTION 

The interpolation of ancillary meteorological and ozone 
data used for level -2 processing is made in two steps. Grid- 
ded ancillary counts are spatially and temporally approx- 
imated for all pixels within a SeaWiFS swath. The inter- 
polation in space uses rectangular bilinear approximation. 
The interpolation in time follows the spatial interpolation 
and applies a simple linear approximation. 

The ozone interpolation employs gridded equidistant 
cylindrical maps of ozone obtained by TOMS and the Tele- 
vision Infrared Observation Satellite (TIROS) Operational 
Vertical Sounder (TOVS) instruments. Previously, all an- 
cillary data files [NCEP, TOMS ozone, and the TOVS 
ozone] were assumed to be generated as global gridded files 
for fixed epochs, so that the entire file could be assumed 
to have the same time tag. For the TOMS data, a single 
file was generated for each GMT day, so the time tag was 
assumed to be 1200 UTC on the day. In recent discussions 
with the TOMS Project, it was learned that the TOMS 
ozone data are actually generated in a similar manner to 
the SeaWiFS level -3 binned products, in that the data are 
gridded for whatever time (or times) the sensor viewed the 
locations on that day. In addition, the TOMS satellite 
orbit has a local noon descending node, like that of Sea- 
WiFS, so the observation times for SeaWiFS and TOMS 
are fairly close (within approximately 1/2 orbit). 

Based on the above, it was decided to change the TOMS 
ozone data interpolation to reflect the actual observation 
times for the gridded data, which varies linearly from 
0000 UTC at 180° longitude, to 0000 UTC at -180°. The 
interpolation for the other ancillary data types remains un- 
changed. Figure 41 shows an example of a TOMS daily 
map with the latitude and longitude, and starting and 
ending times for the four corners of the grid. The lat- 
itudinal distance between consecutive grid points is 1.0° 


and the longitudinal distance between the grid points is 
1.25°. Ozone maps are, therefore, composed of 180 points 
in the vertical direction and 288 points in the horizontal 
direction. 

Three TOMS or TOVS files are normally submitted 
for spatial and temporal approximation of ancillary ozone 
data for a whole SeaWiFS swath. Two files out of the three 
are used at the ancillary ozone approximation of each Sea- 
WiFS scan line. This chapter describes the modification 
introduced to the temporal approximation of TOMS ozone 
ancillary data. The modification calculates the true time 
for each ozone grid point based on information provided 
within a TOMS file header. 

10.2 NEW OZONE SCHEME 

The linear function of time in terms of longitude for 
TOMS ozone files is given in Fig. 42. The function can be 
written as 


T t 


(T s — T e )X T s 4- T e 
360.0 2.0 


(37) 


where Tp is the (sequential day of the year) time of a 
TOMS ozone observation along a given longitude, X is the 
corresponding longitude, and T s and T e are the (sequential 
day of the year) start, and end times, respectively, of the 
TOMS file. 

The modified temporal approximation of TOMS ozone 
ancillary data applies the difference, AT, between the time 
of the SeaWiFS observation of the current scan line and the 
TOMS time associated with the corresponding longitude 
for each scan line pixel (AT^): 


A T n = \T C - TtU (38) 
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Min : -51.4200 Mean : -1.0137 

Max: 55 6800 STD : 9.8157 



Min : -40.6800 Mean : -0.5301 

Max: 20.5000 STD : 7.2945 


Fig. 43. Histograms of ozone differences between the old and new schemes of TOMS approximation for 
SeaWiFS swaths captured around noon (left panel) and the IDL (right panel). 


where T c is the time of SeaWiFS capture of the present 
scan line. The same time, T c , is used for all pixels within 
a scan line, because the fraction of a second which it takes 
SeaWiFS to attain the complete scan line can be neglected 
in the calculations. 

The old temporal approximation defined A T as a dif- 
ference of scan line time from the beginning of the TOMS 
file: 

A T° = ITc-T.I, (39) 

Although the old interpolation scheme (AT 0 ) assumed the 
T s to be 1200 UTC for an entire daily ozone map, it was 
reading the actual T s value from the TOMS file header 
which returned the time around 0000 UTC. 

Temporal approximation computes weights for ozone 
grid points on the corners of the cell which spatially en- 
compasses a SeaWiFS scan line point. The weights are 
calculated for both ancillary TOMS files used in the inter- 
polation, W\ and W 2 for the first and second file, respec- 
tively. Two AT values (ATi and AT 2 ) for each file are in 


the same way converted into weights for the old and 

new 

processing: 

Wl ATi + A T 2 ’ 

(40) 

and 

ATi 

W 2 1 

A Ti + A T 2 ' 

(41) 

The final interpolation equation is as follows: 



F, = W l F [ + W 2 F 2 , 

(6) 


where iq is the interpolated ozone value and F\ and F 2 
are spatially approximated ozone values within the TOMS 
first and second file, respectively. 

The new approximation puts the highest weight on 
ozone grid cells from the ancillary file that is most con- 
current with the given scan line point. A low weight is 
given to ozone grid cells from the other ancillary file later 


in time. A switch of ancillary files takes place when a scan 
line crosses the International Date Line (IDL). The old ap- 
proximation algorithm allocates a higher weight to ozone 
grid cells from the ancillary file whose start time is the 
closest to the scan line time which gives greater weight to 
the nonconcurrent ozone counts at the start of the TOMS 
map. 

Coincident TOMS and SeaWiFS acquisition times re- 
sults in progressively low AT for ozone grid cells from the 
ancillary file that is most concurrent with scan line points. 
The new algorithm works well for all ancillary file data cap- 
tured over a period of time, such as TOMS ozone where 
T e — T s ~ 24 h. It is also appropriate for ancillary data 
where the common time for the whole file is mid-time, T m , 
and T m = T e = T s . The general concept of the method, 
therefore, could also be applied to other ancillary files, 
such as TOVS. Currently, the new processing only acti- 
vates when there are two ancillary files submitted to the 
approximation of a scan line and both are TOMS ozone. 

10.3 RESULTS 

The new TOMS ozone approximation scheme intro- 
duces noticeable differences for SeaWiFS swaths in both 
ozone content and chlorophyll concentration. The mod- 
ified algorithm was executed on eight-day data and pro- 
duced level -2 products. A two-day run produced level -2 
and quality control products for both conventional and al- 
tered algorithms of TOMS ozone approximation. 

10.3.1 Ozone Comparison 

Old and modified TOMS ozone approximation schemes 
were compared. There was a substantial difference in ozone 
counts both with over- and underestimation. Figure 43 
compares histograms based on the ozone value differences 
between the old and new schemes for a near-noon and near- 
midnight (near the IDL) SeaWiFS swaths. Within the 
entire two-day sequence of SeaWiFS processing, the new 
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Chlorophyll Concentration 
(old ozone approximation algorithm) 



Chlorophyll Concentration 
(old ozone approximation algorithm) 


Fig. 44. Scatterplots of two-day chlorophyll concentration values obtained by the old approximation scheme 
as a function of the chlorophyll level difference between the old and new algorithms. 


method of ozone approximation from TOMS files produces 
a difference of up to 20% from the previous algorithm. 

The switch of the concurrent TOMS files along the IDL 
results in discontinuities in the ozone approximation val- 
ues. The artifacts are created because successive day-to- 
day ozone maps do not continuously fall one into another, 
and the new scheme puts the largest weight on ozone cells 
located in the most concurrent file causing the switch at 
the IDL. The amount of the discontinuity reaches up to 
12-17%. 

10.3.2 Chlorophyll Comparison 

The modified algorithm of ozone approximation from 
TOMS files has caused measurable changes in the chlo- 
rophyll concentration values compared to the old method. 
The results were obtained on SeaWiFS GAC data captured 
over two days and processed using ozone information cal- 
culated according to the old and new scheme of TOMS 
interpolation. In the comparison of old and new chloro- 
phyll concentration results, only those pixels which were 
not flagged by either level-2 processing algorithm were con- 
sidered. 

Figure 44 shows two scatterplots of differences in chlo- 
rophyll levels produced by the old and new methods as a 
function of chlorophyll concentration using the old ozone 
approximation scheme. The first plot shows the full scale 
of chlorophyll levels and data differences, and the second 
plot is a magnification of the smaller chlorophyll concen- 
trations and difference values. The scatterplots indicate 
that there are some isolated pixels for which the difference 
in chlorophyll levels between the two methods is as high 
as 100%. 

Histograms of chlorophyll concentration differences be- 
tween the two schemes for a two-day period of SeaWiFS 


capture are displayed in Fig. 45. The same data are shown 
for two ranges of differences in chlorophyll levels. The per- 
centage of significantly changed chlorophyll values is small, 
therefore, only narrow domains around the zero level of 
chlorophyll difference can show on the histogram. 

Overestimation of chlorophyll levels by the old algo- 
rithm is much more common, about 72%, for larger chlo- 
rophyll difference values exceeding 0.1 mg m -3 . Within the 
lower difference range (below 0.1 mg m" 3 ), 74% of counts 
are underestimated by the old method in comparison with 
the new scheme. 

The overwhelming majority (96.8%) of chlorophyll con- 
centration variances between the two methods account for 
chlorophyll level differences smaller than 0.1 mgm“ 3 . This 
value could be nonetheless significant, as the chlorophyll 
concentration of much of the global ocean is in the range 
of 0.1 mg m 3 . The bulk of chlorophyll changes, 70%, is 
below 0.01 mgm~ 3 . 

The discontinuity of new ozone approximation along 
the IDL is not visually discernible even in rescaled chlo- 
rophyll concentration images. The difference can only be 
observed in comparison with chlorophyll counts obtained 
with the old method of ozone approximation. Within the 
limited data set, only intermediate chlorophyll concentra- 
tion levels around 0.3 mg m -3 were observed for the IDL 
regions. The differences between the two algorithms in 
these areas can be estimated up to 20% for either side of 
the IDL. 

10.4 CONCLUSIONS 

A new algorithm for the temporal interpolation of ancil- 
lary TOMS ozone data was needed to reflect the actual ob- 
servation times for TOMS gridded maps. The new scheme 
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Data Values of Chlorophyll [mg nr 3 ] Data Values of Chlorophyll [mg nv 3 ] 

Min : -39.2850 Mean : -0.00143761 

Max : 56.4080 STD : 0.244098 


Fig. 45. Histograms of differences in chlorophyll concentrations between chlorophyll levels obtained using 
the old and new ozone approximation schemes for two days of SeaWiFS data. 

introduces measurable differences in chlorophyll concen- the change is below 0.01 mg m -3 . Within the 3.2% range 
tration counts compared to the old algorithm. In 96.8% of cases, there are some isolated occurrences when the new 
of the cases considered, however, the change in chlorophyll method causes a large variation in chlorophyll readings 
concentration is below 0.1 mg m~ 3 , and in 70% of the cases, reaching as much as 100%. 
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CVT Calibration and Validation Team 
CZCS Coastal Zone Color Scanner 
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HPLC High Performance Liquid Chromatography 
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IDL Interactive Data Language or International Date 
Line (depending on usage). 

JRC Joint Research Centre (Ispra, Italy) 

LAC Local Area Coverage 

MOBY Marine Optical Buoy 

MOCE Marine Optical Characterization Experiment 
MODIS Moderate Resolution Imaging Spectroradiome- 
ter 

MOS Marine Optical Spectroradiometer 

NCEP National Center for Environmental Prediction 
NCSA National Center for Supercomputing Applica- 
tions 

NDVI Normalized Difference Vegetation Index 
NIR Near-Infrared 

OCTS Ocean Color and Temperature Sensor (Japan) 

PlyMBODy Plymouth Marine Bio-Optical Data Buoy 
PROSOPE Productivity of Pelagic Oceanic Systems (French 
translation) 

QC Quality Control 

RH Relative Humidity 


SeaBAM SeaWiFS Bio-optical Algorithm Mini- workshop 
SeaBASS SeaWiFS Bio-optical Archive Storage System 
SeaDAS SeaWiFS Data Analysis System 
SeaPRISM SeaWiFS Photometer Revision for Incident Sur- 
face Measurement 

SeaWiFS Sea- viewing Wide Field-of-view Sensor 
SIMBIOS Sensor Intercomparison and Merger for Biologi- 
cal and Interdisciplinary Ocean Studies 

SIO Scripps Institution of Oceanography 
SIRREX SeaWiFS Intercalibration Round-Robin Experi- 
ment 

SQM SeaWiFS Quality Monitor 
SXR SeaWiFS Transfer Radiometer 

TAO Tropical Atmosphere-Ocean 
TIROS Television Infrared Observation Satellite 
TOA Top of the Atmosphere 
TOGA Tropical Ocean Global Atmosphere 
TOMS Total Ozone Mapping Spectrometer 
TOVS TIROS Operational Vertical Sounder 

UTC Coordinated Universal Time (definition reflects 
actual usage instead of following the letters of 
the acronym) 

Symbols 

a 0 Constant term of function of illuminated frac- 
tion of the lunar surface. 

a\ Linear term of function of illuminated fraction 
of the lunar surface. 

6 0 Constant term of interpolation of lunar phase 
function. 

b\ Linear term of interpolation of lunar phase func- 
tion. 

b 2 Quadratic term of interpolation of lunar phase 
function. 

Ci (A) Wavelength-dependent phase angle correction 
factors. 

C a Chlorophyll a. 

Cdark Dark count from sensor output data. 

Gout Counts from sensor output data. 

d Detector. 

Dim SeaWiFS-moon distance in the mean radius of 
the lunar orbit (R M ). 

Dsm Sun-moon distance in astronomical units (AU). 

E s Surface irradiance. 

/i(A) Illuminated fraction of the lunar surface. 

/2(A) Interpolated lunar phase function. 

/3(A) Difference between measured and fitted inte- 
grated lunar radiances. 

/4(A) Correction factor for the variation in the az- 
imuthal angle. 

Fq(X) Extraterrestrial solar irradiance adjusted for the 
Earth-sun distance variations. 

F\ Spatially approximated ozone count for the first 
ancillary TOMS file. 

F 2 Spatially approximated ozone count for the sec- 
ond ancillary TOMS file. 

Ft Final interpolated value of approximated ozone 
count. 

f /Q Ocean surface bidirectional reflectance parame- 
ter. 
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g Gain. 

Gr Gain factor for band 7. 

Gs Gain factor for band 8. 

£ 7:8 Combined gain for band 7 and 8. 

h Shorthand notation variable, see (19). 

K Approximated aerosol effects in atmospheric dif- 
fuse transmittance. 

K\ Counts-to-radiance conversion factor. 

K 2 Detector temperature-dependent correction fac- 
tor. 

K 3 Scan modulation correction. 

N(490) Diffuse attenuation coefficient at 490 nm. 

K a Aerosol effects (0.0054) estimated for an aerosol 
optical thickness of 0.1 and the maritime aerosol 
model for 90% relative humidity. 

L( A) Radiance. 

L a (A) Radiance measured at the TO A from aerosol 
scattering alone. 

La Aerosol radiance equal to L a + L ra . 

Lb Lwn ratio constructed from band A divided by 
band B. 

L os ( A) Radiance arising from aerosol single scattering. 

L C (A) Fitted integrated lunar radiance. 

L/(A) The foam, or whitecap-leaving, radiance at the 
sea surface. 

L g ( A) Sun glint radiance at the sea surface. 

Lgn(A) The normalized sun glint radiance. 

L m (A) Measured integrated lunar radiance. 

L r (A) Radiance measured at the TOA from Rayleigh 
scattering alone. 

L ro (A) Radiance at the TOA from Rayleigh- aerosol in- 
teractive scattering. 

Ls(A) Calibrated at-sensor radiance. 

Lt Sun glint corrected radiance. 

Lt(A) Radiance measured at the TOA. 

Lw(A) Water-leaving radiance. 

Lvvn (A) Normalized water-leaving radiance. 

Lwn The ratios of the SeaWiFS values to the MOBY 
values for normalized water-leaving radiance. 

M Half-angle mirror side correction factor. 

ms Half-angle mirror side. 

Nm Mean number of scan lines in a lunar image. 

N\ Lunar radiance normalizing factor to a common 
sun-moon distance. 

N 2 Lunar radiance normalizing factor to a common 
SeaWiFS-moon distance. 

N 3 Lunar radiance normalizing factor for the illu- 
minated fraction of the lunar surface. 

N 4 Lunar radiance normalizing factor to a common 
number of scan lines in a lunar image. 

N 5 Lunar radiance normalizing factor to a common 
phase angle. 

Ne Lunar radiance normalizing factor for the wave- 
length-dependent phase angle correction. 

pxl Pixel number along a SeaWiFS scan line. 

r Ratio between two aerosol models from the de- 
rived 6 values. 

R Distance between two bodies. 

R b Ratio of the SeaWiFS diffuser BRDF to lamber- 
tian BRDF. 

R m Mean radius of the lunar orbit. 


t Time tag of the sensor output data. 
f(A) Diffuse transmittance of the atmosphere at the 
sensor viewing direction. 

to Reference time for temporal corrections to the 
T Detector temperature from the sensor output 
data, radiometric response of the instrument. 
T(A) Direct transmittance of the atmosphere at the 
sensor viewing direction. 

To(A) Direct transmittance of the atmosphere at the 
solar direction. 

T c Time of SeaWiFS capture of a scan line. 

T e Sequential day of the year time for the end of 
TOMS file. 

Tm Mid-time. 

Tref Detector reference temperature. 

T 3 Sequential day of the year time for the start of 
TOMS file. 

T t Sequential day of the year time of TOMS ozone 
observation along a given longitude. 

W\ Weights used in the temporal approximation of 
ozone counts for the first ancillary TOMS file. 

W 2 Weights used in the temporal approximation of 
ozone counts for the second ancillary TOMS file. 

X Longitude. 

a SeaWiFS vicarious gain. 

c*(A», Aj) Angstrom coefficient for the wavelengths A, and 

Aj* 

a°(At, 865) Angstrom exponent from the derived model a. 
a 6 (Aj, 865) Angstrom exponent from the derived model b. 

p Constant term in the temporal correction to the 
radiomteric response of the instrument. 

7 Linear term in the temporal correction to the 
radiometric response of the instrument. 

6 Quadratic term in the temporal correction to the 
radiometric response of the instrument. 

AT Difference of scan line time from the beginning 
of the TOMS file. 

ATi AT value converted into weights for the old and 
new processing. 

AT 2 AT value converted into weights for the old and 
new processing. 

A T n Difference in Julian time between SeaWiFS and 
TOMS observations for the new algorithm. 

A T° Difference in Julian time between SeaWiFS and 
TOMS observations for the old algorithm. 

£(765,865) Atmospheric correction parameter. 

£(A, 865) £ value of atmospheric correction for wavelengths 
A and 865 nm. 

€ (Ai, A j) Ratio of single-scattered aerosol reflectance be- 
tween two wavelengths, p as (Ai)/pas(Aj). 

0 Sensor zenith angle. 

0o Solar zenith angle. 

1 ? Phase angle of the lunar calibrations. 

A Wavelength. 

A*,Aj Wavelengths corresponding to sensor spectral 
bands i and j. 
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p Reflectance. 

pa (A) Reflectance measured at the TOA from aerosol 
scattering alone. 
pA Aerosol reflectance. 
po*(A) Single-scattered aerosol reflectance. 

P/( A) Reflectance at the sea surface from foam (or 
whitecaps). 

p r (A) Reflectance measured at the TOA from Rayleigh 
scattering alone. 

Pro (A) Reflectance at the TOA from Rayleigh- aerosol 
interactive scattering. 
pt(A) Reflectance measured at the TOA. 
pw{ A) Water-leaving reflectance at the sea surface. 

a Standard deviation. 

<t(X) Standard deviation of quantity X. 

r 0 (A) Aerosol optical thickness. 

To (At) Initial estimate for r a (A). 
t 0 *(A) Ozone optical thickness. 
r r (A) Rayleigh (air molecules) optical thickness. 

<t> Azimuthal angle. 

<p Azimuthal angle of incident sunlight on the Sea- 
WiFS diffuser 

Angle of rotation in azimuth relative to the Sea- 
WiFS diffuser. 
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